SEMICLASSICAL LOW ENERGY SCATTERING FOR ONE DIMENSIONAL 
SCHRODINGER OPERATORS WITH EXPONENTIALLY DECAYING 

POTENTIALS 



OVIDIU COSTIN, ROLAND DONNINGER, WILHELM SCHLAG, AND SALEH TANVEER 



Abstract. We consider semiclassical Schrodinger operators on the real line of the form 



with h > small. The potential V is assumed to be smooth, positive and exponentially decaying 
towards infinity. We establish semiclassical global representations of Jost solutions f±{-,E;h) with 
error terms that are uniformly controlled for small E and h, and construct the scattering matrix 
as well as the semiclassical spectral measure associated to H{h). This is crucial in order to obtain 
decay bounds for the corresponding wave and Schrodinger flows. As an application we consider the 
wave equation on a Schwarzschild background for large angular momenta £ where the role of the 
small parameter h is played by It follows from the results in this paper and 0, that the decay 
bounds obtained in [8], [6] for individual angular momenta £ can be summed to yield the sharp t~'^ 
decay for data without symmetry assumptions. 



1. Introduction 

In this paper we study Schrodinger operators H{h) of the form 

Hin)fix) := -h^f'ix) + V{x; h)f{x) (1.1) 

for X G M. The potential V is assumed to be positive, smooth and exponentially decaying as 
|x| — )• cxD. As indicated by our notation, we also allow V to depend on the semiclassical parameter 
but only in a very mild way, i.e., all assumptions on V are assumed to hold uniformly in small h. 
This will be important for our main application, the wave equation on a Schwarzschild background. 
More precisely, we require V{-;h) G C°°(M) to be of the form V{x;h) = V±{e~^^^;h) for itx > a 
where a > is some (large) constant and V±{-]fi) are smooth functions with V±(0; K) = 0. It follows 
from the standard theory (see e.g., [20], [21]) that under the above assumptions the operator H{fi) 
has a self-adjoint realization on L^(M) with domain T){H{h)) = ff^(M) and its spectrum is purely 
absolutely continuous with a{H{K)) = [0, oo). The present paper is devoted to the study of low 
energy {E — )• 0+) scattering in the semiclassical limit h — )• 0+, i.e., we consider the low energy 
and semiclassical limits simultaneously. To be more precise, we are primarily interested in global 
semiclassical representations of the Jost solutions f±{x,E; K) and their derivatives with respect to 
E that hold uniformly in small E and h. As a consequence, we obtain uniform bounds on relevant 
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quantities that are constructed based on the Jost solutions such as the reflection and transmission 
amplitudes as well as the spectral measure. 

Our main motivation for investigating this problem comes from the large angular momentum be- 
havior of solutions to the Regge- Wheeler equation which describes perturbations of a Schwarzschild 
black hole in general relativity [15]. The Regge- Wheeler potential decays exponentially towards 
the event horizon of the black hole and thus, one is naturally led to a Schrodinger operator of 
the form (jl.ip where corresponds to the angular momentum i. Thus, in the limit £ — )• (X) 
the problem becomes semiclassical. Consequently, the present paper is part of the series [8], [7], 
[U] devoted to the study of wave evolution on a Schwarzschild background. However, we believe 
that our results and techniques are of independent interest and therefore we decided to embed our 
work in the broader context of rigorous quantum mechanical scattering theory in the semiclassical 
limit. Scattering theory is by now a classical subject and we cannot possibly do justice to the 
vast literature in that area. We confine our discussion of the existing literature to selected works 
devoted to the one-dimensional case which are related to our problem. 

A standard reference for scattering on the real line is [6] which provides many of the by now 
classical results on the subject. Low energy scattering for one-dimensional Schrodinger operators 
is studied in, e.g., [10], [12], [2\, [22j whereas the semiclassical limit is considered in However, 
we find that the available results are not suitable for our purposes. The only paper we are aware 
of that studies the double asymptotics E,h ^ 0+ in a way similar to us is which deals with 
inverse square potentials. Since the Regge- Wheeler potential exhibits inverse square decay towards 
spatial infinity, [4J is used to deal with the far field region in our main application which we now 
describe in more detail. 

1.1. The problem of wave evolution on a Schwarzschild background. One of the major 
open problems in mathematical general relativity is a proof of the stability of black holes as de- 
scribed by the Kerr solution. As a first step in this direction one considers the behavior of linearized 
perturbations of Schwarzschild black holes. It is well-known that such perturbations are described 
by the Regge- Wheeler equation 



which is effectively a one-dimensional wave equation with a potential. The number £ G No denotes 
the angular momentum and the coordinate system (t, x) is chosen in such a way that x — >• oo 
corresponds to spatial infinity whereas x — )• — oo is the location of the event horizon of the black 
hole. Furthermore, the Regge-Wheeler potential Ve reads 



where a E {—3,0, 1} is a fixed parameter which corresponds to different physical situations. The 
function r(x) is implicitly defined by the equation x = r{x) + log(r(2;) — 1). Consequently, Vi{x) 
behaves like ^^^i^ as x — )■ oo and decays exponentially as x ^ — oo. More precisely, for x < 0, Ve{x) 
can be written as a convergent series in powers of e^, see [7j. We refer the reader not familiar with 
this problem to the standard literature, e.g., [3] or the introduction of [8j as well as the references 
therein. Recently |8] , [6] established decay bounds for solutions of Eq. (|1.2p of the form 



\\{r'Mtr)\\L^<C,{t)-^\\{-)'MOr)\\L^ + \\{-)'^^^^ (1.3) 



for t > and any ^ G Nq. Since the Regge-Wheeler equation ()1.2p results from an angular 
momentum decomposition, the ultimate goal is to sum the decay bounds p.Sp over all i. However, 
in order to conclude the overall decay one needs to control the behavior of the constants Cg 
in (II. 3p as £ — 7- oo. At this point the semiclassical analysis of the present paper (and also [4J) 
comes into play. Indeed, by dividing Eq. (II. 2p by h~'^ := £{i + 1) one is led to a semiclassical 



d^Mt,x)-dl 



ipeit, x) + V(>{x)^e{t, x) = 0, {t, x) € [0, oo) x M 



(1.2) 
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problem of the type (jl.ip . As a consequence, the Regge- Wheeler potential is a prominent example 
from mathematical physics where our results apply, at least for x < 0. In the case x > 0, where 
the potential decays like an inverse square, one has to rely on [4j. In fact, in \7] it is shown how 
to synthesize [3] and the present results to obtain the sharp decay for linear perturbations of 
Schwarzschild without symmetry assumptions on the data. At this point it has to be remarked that 
the long-sought decay also follows from an independent result by Tataru [T8] which appeared 
simultaneously to [7J. An important consequence of our results (in conjunction with [4j) is that 
the small energy contributions decay rapidly as £ — )• oo and do not present an obstruction to the 
summation. Indeed, as expected, the important contributions come from energies close to the 
maximum of the potential [7J. 

1.2. Notations and conventions. Throughout this work we write (•) : M — t- [^,oo) to denote a 
smooth, symmetric function that satisfies (x) = \x\ for > 1. Furthermore, we assume (•) to be 
strictly monotonically increasing on (0, oo). For two real numbers a, b we write a < 6 if there exists 
a constant c > such that a < cb. Unless stated otherwise, it is assumed that c is absolute, i.e., 
it does not depend on any of the quantities involved in the inequality. Similarly, we use a>b and 
a c^b means a <b and 6 < a. Asymptotic equality is denoted by the symbol ~. For the Wronskian 
W{f,g) of two functions f,g we use the convention W{f,g) := fg' — f'g. The expression 0{f{x)) 
for some function / is used to denote a generic real-valued function that satisfies |0(/(x))| < |/(a;)| 
in a domain of x that follows from the context. We write Oc{f{x)) if the functions attains complex 
values as well. Finally, the letter C (possibly with indices) stands for a positive constant that may 
change its value at each occurrence. 

1.3. The scattering matrix. The (semiclassical) Jost solutions f±{x, E; K) of the operator H{h), 
where E > 0, are defined by the relation 

H{h)f±{;E;h)=E^f±{;E;h) 

and the asymptotic conditions f±{x,E;H) ~ e^*"ft^ as x — )■ ±00. Note, however, that f±{-,E;h) 
are not eigenfunctions of H[h) since they do not belong to L'^(R). It is a classical result (see 
e.g., [5]) that the Jost solutions exist if V{-;h) € L^(M) (which is clearly the case here) and that 
they are unique. The physical significance of the Jost solutions emerges from the fact that they 
correspond to outgoing waves. This terminology is explained by the observation that the functions 

ip±{t, X, E; h) = e~'^~^ f±{x, E;h) satisfy the Schrodinger equation 

ihdt^±{t,-,E;h) = H{h)iP±{t,-,E;h) 

and they behave asymptotically as 'ip±{t,x, E;h) ~ e-*7r(^*T^) for x — t- ±00. In other words, 
they describe waves that travel towards ±00. It should again be remarked here that traditionally 
one writes E instead of E^ but in order to avoid square roots we prefer to use E'^. From the 
construction of the Jost solutions via Volterra iterations it also follows that the asymptotics can 
be formally differentiated, i.e., f'^{x, E; h) ~ zti-^e^*"^ as x — )• ±00 and we immediately obtain 
the (x-independent) Wronskians 

W{f±{;E;n),f±{;E;h))=T2i^ 

by evaluation as x — )• ±00. Consequently, if / 0, {f±{-, E; h), f±{-, E; h)} are two fundamen- 
tal systems for the ordinary differential equation H{h)f = E'^f and there must exist connection 
coefficients a±{E;fi) and b±{E;h) such that 

f±{;E; h) = a^{E- h)f^{; E; h) + b^{E- h)f^{;E;h). 

^Note carefully our somewhat unusual convention. We denote by E the square root of the energy rather than the 
energy itself. This turns out to be more convenient. 
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This yields 

where we omit the arguments E and h occasionally in order to improve readability. We obtain 
\b+{E;h)\ > 1 and consequently, 

\W{f-,U)\ = \W{a+U + b+U,U)\ = 2|6+|f > 2^ (1.4) 

which shows that {f^{-, E; h), f+{-, E; h)} and {/-{■, E;h), f^{-, E;h)} are also two fundamental 
systems for the equation H{K)f = E'^ f provided that E ^ 0. This is already a nontrivial statement 
since it contains global information. For obvious reasons these two systems are called outgoing and 
incoming, respectively. If E = 0, {/_(•, 0; h), /+(•, 0; h)} may or may not be a fundamental system 
for H{h)f = 0, depending on the special form of the potential. In the latter case one speaks of the 
existence of a zero energy resonance. However, we are not faced with this complication since it is 
not hard to see that the positivity assumption on the potential already excludes the existence of a 
zero energy resonance, see @]. A general solution f{-,E; h) of H{h)f{-,E; h) = E'^f{-,E; K) can be 
expanded as 

/(x, E; h) = c7 {E; h)f^ (x, E- K) + c+{E- h)U (x, E; h) 

or, alternatively, 

/(x, E- h) = c+{E- h)U{x, E- h) + Co{E; h)f_{x, E; h) 
and there exists a linear transformation that relates the representation with respect to the incoming 
basis (/_, /+) to the one with respect to the outgoing basis (/+, /_). The 2x2 matrix that represents 
this transformation is denoted by §(-E; H) and called the scattering matrix. We write 

V S21 S22 / 

and in the special case cj{E; h) = 1, c'^{E; h) = we obtain 

f-{x,E;h) = sn{E; h)f+{x, E; h) + S2i{E- h)f.{x, E; h). 

The physical interpretation of this relation is that an incoming wave from the left scatters at the 
potential and gets decomposed into a transmitted and a reflected part. Consequently, sn and S21 are 
called transmission and reflection amplitudes and we write sii{E] h) = t{E; h), S2i{E; h) = r{E; h). 
A similar statement is true for an incoming wave from the right. It is not hard to see that S22 = t 
and |ip + |rp = 1 which complies with the probabilistic interpretation of the quantum mechanical 
wave function. Furthermore, it is well-known that S{E; K) is unitary and completely determined 
by t{E] K) and r{E; h). In semiclassical low energy scattering one is interested in the behavior of 
S(i?; K) as both E and h tend to 0, simultaneously and independently of each other. 

Let us remark that in the pure quantum mechanical setting, where /i > is fixed, the behavior of 
E:{E; H) for E — t- 0+ is trivial. This follows from the fact that for exponentially decaying potentials 
the Jost solutions f±{x,E; h) are smooth in E around |5J. On the other hand, if > is fixed 
and only h — t- 0+ is considered then the standard WKB method works and the corresponding 
result is well-known and classical. However, the problem becomes highly nontrivial if one allows 
both limits E,h ^ 0+ and moreover, if one is interested in uniform bounds for all derivatives 
with respect to E. In order to approach this problem, we rely on a coordinate transformation 
which maps H{fi)f = E"^ f to a perturbed Bessel equation. The resulting equation is solved by a 
perturbative construction around Bessel functions which is based on suitable Volterra iterations. 
At this point the asymptotic theory of Bessel functions becomes crucial. For the convenience of the 
reader we have compiled the necessary results on Airy and Bessel functions in the two Appendices 
O and |DJ However, since we need to control all derivatives of the perturbative solutions, the 
standard results on Volterra equations are insufficient and we have to extend them. This is done 
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in Appendices |B] and O where the required theory is developed. As a consequence, we obtain the 
following representation of the Jost function of the operator H{h). 

Theorem 1.1. Fix small constants Eq^Hq > and let H{h) he the semiclassical Schrodinger 
operator on L^(M) given by 

H{h)f{x) := -h^f"{x) + V{x; h)f{x), x G M, /i G (0, ^o) 

where the potential V satisfies 

(1) V^''^ G C(M X [0, ho]) for all G No, 

(2) V{x; h)>0 for allxGR, h£ [0, ^o); 

(3) V{x;h) = V±{e~^^^;h) for itx > a where a > is some fixed constant and the functions 
V± have the property that vl^^ G C([0, 1] x [0, ho]) for all G No and V±{0; h) = for all 
h G (0,/io). 



Then, with a := y^+4:E'^, the outgoing semiclassical Jost function f+(0, E;h) of the operator 
H{h) has the representation 

f+{0,E; h) = aij+{E; /j)ei(^+(^;'^)+*^+(^;'"0) h + nai{E; h) + a2{E; h)e-i^+^^'^^ 
f+{0,E;h) = a2h-^c+{E;h)j+{E;h)erA^+^^'^^+''^+(^'^^') U + ha-i{E;h) + ai{E;h)e-^^+^^'^^ 

where c+{E;h),S+{E;h) > 1, T+{E;h) > E\log{^+4:E'^)\, \'y+{E;h)\ ~ 1, and we have the bounds 

\d^EA{E;h)\ <Cih-\ A G {c+,7+,5+,r+,ai,a2,a3,a4} 

for all I G No, E G (0, Eq) and h G (0, ho). Furthermore, cri{E; h) and ct^{E] h) are real-valued. 

Based on the representation in Theorem II. II we obtain our main result on the scattering matrix. 

Theorem 1.2. Under the assumptions of Theorem M.li the reflection and transmission amplitudes 
r{E; h) and t{E; h) associated to H{h) are of the form 

t{E- h) = t{E; ;,)-^e-^(^(^^'"')+^^(^^'"'))[l + ha,{E; h)] 

r{E- h) = p{E- ;i)e-f + nu2{E- h)] 



for functions S, T, r_, r and p satisfying S{E;h) > 1, T{E;h),T^{E;h) > log(^ + ^E'^)\, 
\t{E; h)\ ~ \piE; ~ 1 for E G (0, £^0) ^.''^d h G (0, Hq). Furthermore, we have the bounds 

Id'^AiE-h)] < Cih-\ A G {S,T,r_,T,p,ai,a2} 

for all I G No, E G (0, Eo) and h G (0, ho) where Ci > are constants that only depend on i. 

Some remarks are in order: 

• Theorem 11.21 shows that the probability for a quantum particle to tunnel through the 
potential barrier decays exponentially as h ^ 0+, i.e., as one approaches the classical 
regime. As already mentioned, for fixed -E > this exponential behavior is well-known, see 

e.g., m- 

• The factor in t{E; h) reflects that zero energy is not a resonance, which is equivalent 
to the vanishing of the transmission coefficient at zero energy. 



'By V^^*^ we mean the fc-th derivative with respect to the first variable, i.e., V^''\x\ h) = d^V{x\ h). 
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Theorem 11.21 yields control over all derivatives with respect to E of the involved quantities. 
This is the most salient feature of the result and it unveils the main difference to Schrodinger 
operators with potentials that exhibit power law decay where one loses powers of E upon 
differentiating with respect to E, see e.g., |4J. Such behavior is in stark contrast to the 
situation here where one only loses powers of h and such a loss is in fact negligible compared 
to the size of e~ s . In particular this shows that the reflection and transmission amplitudes 
are smooth in E at E = which is not the case for potentials that decay according to a 
power law. 

The functions S{E;h), T{E;h) and T-{E;h) are not uniquely determined and we do not 
state explicit formulas here. However, comparison with the classical WKB result for E > 
small but fixed shows that in this case one may take 



SiE; h) = / yJ\V{x;h)-E'^\dx 

where ^xt^±{fi) > are the two solutions of V{x; h) — E'^ = which are bounded i( E > 
is fixed. Similar statements apply to T{E;h) and T-{E;h), cf., e.g., [H] . 

• Note further that the exponential decay of the transmission amplitude t{E; h) for h — )• 0+ 
should not be confused with the exponential decay found in [22]. The result in [22] is valid 
for fixed /i > in the limit E — )• 0+ and has a completely different origin, namely the slow 
decay of the potential considered there. 

• The positivity assumption on the potential can be dropped and replaced by some weaker 
condition. However, one certainly has to exclude the existence of a zero energy resonance 
for the result to hold. 

1.4. The spectral measure. The solution of the Schrodinger equation 

ihdti^{t, ■) = H{K)ij{t, •) (1.5) 

is of course given by ijj{t, ■) = e~ a*^^'*-*^(0, •) where the exponential has to be interpreted according 
to the functional calculus for self-adjoint operators. It is a consequence of Stone's formula (see 
[20]) that the kernel e~h^^'^^\x,x') of the solution operator to the Schrodinger equation (jl.Sp is 
given by the oscillatory integral 

e-^nmh\x,x') = -4^ / e-Ti'''\{x,x',E;h)EdE, x' < x 
vr/i^ Jo 

with the semiclassical spectral measure 

f.{x',E;h)U{x,E;h) 



e{x, x' , E; K) := Im 



:i.6) 



w{f-{-,E-h),u{;E-h)) 

Evidently, it is important to obtain bounds for the derivatives d^^e{x,x' , E; h), I S No, since these 
immediately translate into decay estimates for the time evolution. In this paper we prove the 
following result. 

Theorem 1.3. Under the assumptions of Theorem the semiclassical spectral measure associ- 
ated to the Schrodinger operator H{h) is of the form 

e(0, 0, E; h) = j{E; h)he-^^^^'^'^ 

where 7, S are real-valued functions and S{E; h) >1 for all E G (0, Eq), h £ (0, Hq). Furthermore, 
we have the bounds 

Id^EliE; n)\ < Cfh-', \d'ES{E; h)\ < Ceh-' 
for all E G (0, Eq), h G (0, ho) and £ G Nq. 
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Again, we make several remarks. 

• In some sense the behavior of e(x, x', E;h) at x = x' = is the most important one which 
is why we only consider e(0, 0, E; h) here. Based on our proof it is straightforward to obtain 
corresponding results for the general case x, x' £ M. 

• We emphasize again that derivatives with respect to E cost powers of h instead of powers 
of E as is the case for potentials that decay according to a power law. This behavior 
translates via stationary phase arguments (cf. [TB], [Ej) into a decay of the corresponding 
time evolution that is faster than for any iV G N. 

• Theorem 11.31 shows that the spectral measure decays exponentially in the semiclassical limit 
h — >■ 0+. Furthermore, the same is true for all S-derivatives of e{0,0,E; h). Consequently, 
the small energy contributions to bounds for the associated Schrodinger flow decay rapidly 
as h ^ 0+. 

2. Preliminary transformations 



As outlined in the introduction, we intend to study the Jost function /+(•, E; K) which is uniquely 
defined by H{h)j\{-,E]h) = E'^ f^{-, E; h) and the asymptotic condition f^{x,E;h) ~ e*"ft^ as 
X — )• oo. By applying the variation of constants formula, it immediately follows that fj^{-,E]h) 
satisfies the Volterra equation 

t-oo 

f^{x,E;h)=e'T- + ^^ / sm{^^{y-x))V[y;h)f+[y,E;h)dy. 

J X 

Evidently, this representation does not give us any useful information in the limit h — t- 0+ which 
we are interested in. This is just an instance of the fact that the potential V{x; h) for large |x| 
cannot be treated as a perturbation in the semiclassical limit h — t- 0+. 

2.1. The Liouville— Green transform. The main tool we are going to use is the Liouville-Green 
transform, see e.g., [13j. For the sake of completeness, we briefly review the main ideas. Consider 
a semiclassical problem of the form 

-h^f"{x) + Qix)f{x) = (2.1) 

for X € / C M, / an open interval. Here, the function Q is assumed to be sufficiently smooth. 
Now suppose 99:/— T-JcMisa diffeomorphism (not necessarily orientation-preserving) onto 
some open interval J and set g{ip{x)) := \ip' {x)\^ f (x) . A straightforward computation shows that 
g satisfies the equation 

^2 /// / Qix) , , ^^ ,2/^3v3"(x)2 \ in'" (x)\ , , , , 



if and only if / satisfies Eq. (12. ip . Observe that, if behaves reasonably well on I, this transfor- 
mation introduces a right-hand side that is of size f? globally and may therefore be treated as a 
perturbation. The point where this procedure turns into an art form is when it comes to choose 
the diffeomorphism ip. This has to be done on a case by case basis. For instance, if Q does not 
vanish on /, one may choose ip in such a way that = sign Q{x) which leads to the classical 

WKB ansatz. It is a historical curiosity that the WKB method is in fact a special case of the much 
older Liouville-Green transform. Another situation which we are going to encounter is when Q 
has a single nondegenerate zero, i.e., Q{xt) = and Q'{xt) 7^ for some xt G /. In this case, xt 
is called a turning point of the equation since the behavior of solutions changes from oscillatory 
to exponentially increasing / decreasing. In order to capture this transition, one chooses 99 such 
that ^i^'i = ^{x). This leads to the so-called Langer transform [llj where the problem at hand 
is mapped to a perturbed Airy equation (see Appendix [U]). However, we remark that the "right" 

7 



choice of (p is not always obvious since it may be advantageous to move some terms of Q to the 
right-hand side of the equation prior to applying the Liouville-Green transform. 



2.2. Reduction to a perturbed Bessel equation. The semiclassical problem we study is given 
by 



h^f"{x) + [V{x- h) - E^]f{x) = 



(2.3) 



for X > xq where G K is fixed throughout. According to the requirements in Theorem 1 1.2 1 we may 
assume that V{x;h) = e~^[l + e{e~^;h)] for x > xq where \e{e~^;h)\ < and e{-;h) is smooth 
around and satisfies e{-;h) > —1 for all < <C 1 (the last condition follows from the positivity 
assumption on the potential). Since we only consider energies E close to 0, Eq. (|2.3p has a unique 
nondegenerate turning point Xt{E; h) > defined by the relation V{xt{E; h);h) = E"^. Assume for 
the moment that V{x]H) = exactly for x > xq. In this case, the turning point is given by 
xt{E; h) = — 21og£' and thus, it is not smooth in ii^ at £^ = 0. Since we want to control derivatives 
with respect to E, it is desirable to have the dependence of the turning point on E be as regular 
as possible. Consequently, it is advantageous to introduce new variables. Based on formula (j2.2p . 



it is reasonable to look for a diffeomorphism (p satisfying (p'{x)'^ = e ^ and thus, (p{x) = 2e 2 
a possible choice. Obviously, (p maps (xo,oo) to (0, yo) diffeomorphically where yo := 2e 



2 



IS 

and 



f'ix) < for all X £ (xo,oo), i.e., ip> is orientation-reversing. 
Eq. (|2.3p transforms into 



Setting /((/j(x)) := \ip{x)\2 f{x), 



n'f'iy) + 



1 



+ 4E^ 

y2 



fiy) = 



on y G (0, yo) which turns out to be a modified Bessel equation. Indeed, setting f{z) 
yields 



(2.4) 



1 



f{z)=0, v:=2\ 



which is solved by z 1— )• -y/z/±ijy(z), the modified Bessel functions. This motivates the following. 
Lemma 2.1. Let (p{x) := 2e~2 and yo := 2e~ 2 be as above. Then the function f, defined by 



f{(p{x)) = \ip'{x)\2f{x) 



- f{x), satisfies the equation 



n'f"iy) + 



l + e{y^-h) 



4 



4E' 



f{y) = 



(2.5) 



for y £ (0, yo) if and only if f satisfies Eq. 
smooth around and we have the bounds 



\d^e{y''-h)\<Cky 



2.3p for X £ {xo,oo). Here, the function e{-;h) is 



max{2-fc,0} 



for all y G (0, yo), < 7i < 1 and k G Nq. 
Eq. ()2.5p with the asymptotic behavior 



Finally, there exists a unique solution f-^{-,E;h) to 



{y ^ 0+) 



f+iy,E;h)r. 

and the outgoing Jost function of the operator H{H) is given by 

f+{x, E; h) = y{x)\-\h{p{x),E- h) = ef /+(2e-t, i?; h). 
Proof. By applying formula (|2.2p . Eq. (|2.3p transforms into Eq. ()2.5p . Furthermore, as already 



remarked, Eq. ()2.3p has a unique solution /+(•, E; H) (the Jost function) that satisfies /+(x, E; K) ~ 
e*"^^ as X — 00, see [5]. Consequently, by setting f^{x,E; h) = \ip' [x^y^'i fj^{ip{x),E:, h) we see that 
f^{-,E;h) is a solution to Eq. (|2.5p and the stated asymptotics of f^{-,E;h) follow immediately 



from the asymptotics of f^{-,E; h) and the expHcit form of (p. The bounds on e are a consequence 
of the assumed properties of V. □ 

From now on we use the shorthand notation := ^ + 4i?^. We are interested in the behavior 
of Eq. (j2.5p for small a. However, this is a singular perturbation problem and it is thus natural to 
rescale the independent variable by setting g{z) := f{az) which transforms Eq. (j2.5p into 

- hlg"iz) + [1 - ^ + eia^z''; h)] ~g{z) = (2.6) 

for z £ (0, ^) where hi := ^. Note that < a < 1 and < h <. 1 implies hi G (0, 2) and thus, 
hi is not necessarily small. The point of this rescaling is that \e{a'^z'^;h)\ < a^z^ and therefore, 
Eq. (j2.6p is a regular perturbation of a modified Bessel equation, at least for small z. Clearly, z 
can become as large as and in this case, the perturbation e is of order 1. However, as we 
will demonstrate now, the perturbation can be made globally small by a suitable Liouville-Green 
transform. 



3. The normal form reduction 
First, we study a useful diffeomorphism that will appear frequently in the following. 
Lemma 3.1. Let ( be defined by 



({x) := sign(x — 1) 



l-^\du 



x>0. 



Then C : (0, oo) R is an orientation-preserving diffeomorphism that satisfies C(l) = and 

2 1 

~ X3, Ci^) — for X > 2. Furthermore, we have the derivative bounds 

\C^''Hx)\<Ck{x)l-'' 
for all X > -^j k £ Nq and the asymptotic behavior 

i[-C(x)]i : 



log X + log 2 — 1 + X eo{x), 



\el^\x)\<Ck 



for all X G (0, |). 

Proof. It is obvious that ( is smooth on (0, 1) U (l,oo) and smoothness around 1 is proved by a 
Taylor expansion, cf. the proof of Lemma 13.21 below. This also yields C'(^) > for all a; > and 
thus, C : (0, oo) — )• M is injective and orientation-preserving. Obviously, ( is also surjective and 
thus, a diffeomorphism. The stated asymptotic behavior of C{x) for x — )• oo is also easily seen 
directly from the definition and the derivative bounds for x > ^ follow from the theory set forth in 
Appendix O For the asymptotic behavior as x — t- 0+ note that, for x G (0, 1), we have 



l[-C{x)]i=l ^l^-ldu = \ogu+ yjl-v?- log (l + 

= — log X — \J\ — x^ + log ^1 + \/l — a 

= — log X + log 2 — 1 + e(x) 

where e is smooth on (—1, 1) and e(x) = O(x^). □ 

Our goal is to "remove" the perturbation e in Eq. ()2.6p by a Liouville-Green transform. Note 
that Eq. (j2.6p has a unique nondegenerate positive turning point zt{p?\ h) which is implicitly defined 
by 



Since e(0;^) = 0, it follows from the smoothness of e{-]K) and the implicit function theorem that 
Zt{-;H) is smooth around and of the form zt{a^; K) = l+0(a^), uniformly in < ^ ^ 1. According 
to Section [2.11 we are looking for a diffeomorphism a^; h) that satisfies 

for z € (0, ^) where ' always means differentiation with respect to the first variable. Clearly, if 
a = 0, we may choose Lp{z, 0; h) = z, and thus, since e{a'^z'^] K) depends smoothly on a^, we expect 
ip{-,a'^\h) to be a smooth perturbation of the identity. Fortunately, Eq. (j3.ip is separable and 
denoting the new independent variable by w (i.e., w = ipiz^c?; h)), integration yields 



l-^\du= / ^|i_ J^ + e(a2n2;n)|dn (3.2) 

which, by montonicity, defines a one-to-one correspondence between w and z and this implicitly 
yields the desired mapping ip{-,a'^;h). We can immediately read off that ip{zt(a'^;h),a'^;h) = 1, 
i.e., the turning point is mapped to 1, and lim2^o+ ¥^(-2, a^; /i) = 0. We also set 

wo{a'^;h):= lim (p{z,a'^;h) (3.3) 

which is finite for any < a ^ 1 and wo{a'^; /i) — )• oo as a — )• 0+. Thus, Eq. (j3.2p defines a bijective 
mapping ip{-, a'^; h) : (0,a~^yo) — ^ (0, wo(«^; ^))- In order to study the analytical properties of 
(f, it is important to realize that ip{-,a'^;h) is smooth on {0, zt{a'^;h)) U {zt{a'^;h),a~^yo) but at 
z = zt{a^]h) (or, equivalently, w = 1) we run into some difficulties since both sides of Eq. (|3.2p 
are nonsmooth there. However, it may still be the case that the two singularities cancel out and in 
fact, this is exactly what happens. In order to see this, compare the function Q from Lemma l3.ll 
with the left-hand side of Eq. (j3.2p which shows that it is natural to raise both sides of Eq. (j3.2p 
to the power |. More precisely, taking the sign issues into account, we obtain 



l-^\du 



t{z, a^; h) 



sigD.{w — 1) I 

where 

t{z, c?\ K) := sign(z — ztia^; h)) 
and it follows that (p is in fact given by ^{z, a'^;h) = C~^{t{z, o?\ h)). 



(3.4) 



3.1. Analytical properties of r. The function t{z, a^; h) diverges for z — )• 0+ and has an appar- 
ent singularity at the turning point z = zt{a'^; h). Furthermore, it also diverges for z — )• a~^yo and 
a — )• 0+. Consequently, we study the analytical properties of r near those points separately and 
we start with the behavior near the turning point. 

Lemma 3.2. Let 6,aQ > be sufficiently small. Then the function t, defined in ()3.4p . satisfies the 
hounds 

\dl,d\T{z, 0?-M< CuA^ - Zt{o?- ;i)|— {l-fc-£,0} 

/or all z ^ (1 - (5, 1 + S), a G (0,ao), < h <. 1 and k,i £ Nq. 

Proof. First of all, we have already noted that the turning point zt (a^ ; ?i) is a smooth function of 
near and satisfies zt{a'^; h) = 1 + O(a^), uniformly in < ^ ^ 1. Consequently, by choosing 
ao > sufficiently small, we can always accomplish that zt{a^;h) € (^5 f) for all < a < ao and 
< ?i ^ 1 which we shall assume in the following and we consider (u, a, h) G (^, |) x (0, ao) ^ (0; ^o) 
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where > is sufficiently small. Now we study the integrand in Eq. ()3.4p . Since the turning point 
is nondegenerate, we can write 

1 - ^ + e{a^u'^;h) = [u - zt{a^;h)]g{u,a^;h), g(zt(a^;h),a^;h) > 1 

with a suitable function g(u,a'^;h) that is smooth with respect to u and in the regime of 
parameters being considered here. Choose an arbitrary N £ N. Since g does not change sign, the 
function g := ^ is smooth as well and a Taylor expansion yields 

N 

g{u, o?\ h) = \J g{u, h) = aj{a'^; h)[u — zt{a^\ h)Y + RNiu, a^; K) 

j=0 

with a remainder of the form 



Rn{u, K) 



Nl 



d^+^~g{s,a^;h){u- s)^ds. 



The coefficients aj (a ; h) are given by 



aj{a^;h) = ijdig{u,a^;fi)\^^-,^(^a2.j,) 

and thus, they are smooth functions of a^. We emphasize that due to the nondegeneracy of the 
turning point we have ao(a^; ^) ^ 1- Furthermore, the remainder satisfies the derivative bounds 

for k,£ £ Nq provided that k + £ < N . Consequently, for u > zt{a'^] h), we obtain an expansion of 
the form 

N 

1 - ^ + e(a2n2; h)=Y, a, (a^; h)[u - zt{a^; h)y+^ 

j=0 

+ [u- zt{a^]h)]'^RN{u,a^;h) 

and integration yields 

N 



v+i 



+ [z- zt{a^; h)]2R^{z, c?\ h) 



provided that z G |) and z > zt{a'^; K) where R^ satisfies the same bounds as R^. We obtain 



r(z, a^; K) 



^[z- zt{a^-K)]2 

'_j=0 
2. fcM D „2 



N 



E 



aj{a , h) 



[z - zt{o?\ h)Y + Rn{z, o?\ h) 



= [z — zt{a ■,h)]RN{z,a ]h) 

where |5^2i9f iiAr(z, a^; < Ck,i provided that k + £ < N and z is sufficiently close to zt{a'^;h) 
(recall that 00(0^; /i) ^1). An analogous calculation for z < zt{a'^]h) shows that the result is in 
fact valid for all z with \z — zt{a'^;h)\ sufficiently small. However, since zt{a'^]h) — )• 1 uniformly 
in ^ as a — )• 0+, we can always find a 5 > such that the above holds for all z S (1 — (5, 1 + 5), 
a S (0,ao) and < h <^ 1 provided that > is chosen small enough. Since N was arbitrary, 
the claim follows. □ 



Next, we consider the behavior as 2; ^ 0+. 
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Lemma 3.3. Let 5,ao > be sufficiently small. Then the function r given by ()3.4p has the form, 

I [-t{z, ; /i)] 5 = - log z + ei {z, ; h) 

where ei satisfies the bounds 

\di2d^,ei{z,a^;h)\<Ck,i 
for all z e (0, 1 - I), a e (0,ao) and <h<^l. 

Proof. Let 6,ao > be so small that Lemma 13.21 is applicable. Moreover, choose ao > so small 
that zt{a^; /i) G (1 - |, 1 + |) for all a G (0, ao) and < < 1. If u e (0, 1 - |), the integrand in 
the definition of r, Eq. (j3.4p . can be written as 



'-ij -1 - e{a^v?; h) = iy^l - - v?e{a^u'^; h) = \ + g{u, h) 
and, since the expression under the square root stays away from zero, g satisfies 

2>^ 



\dl,dtg{u,a^-M<Ck 

for all n G (0, 1 — |), a e (0, oq), < ^ 1 and A;,£ G Nq. We split the integral in the definition 



of T according to 



f 

Jo 



10 

and obtain 



+ 



/ \l ^ ~ ^ ~ ^{ct'^u'^'i K)du = — log z + g{z, h) 

J z 

for z G (0, 1 — |) where g satisfies the same bounds as g. For the second integral note that by 



2 

definition of r we have 



''^ ^^-l-eia^u^;h)du = |[-t(1 - ia'';h)]l 



and according to Lemma 13.21 we have the bounds 

1^.4-^(1 - ia^; h)]l\ < C,\l - I - zt{a';h)\l-' < C, 
since |1 — | — zt{a^; h)\^ 5 for all a G (0, oq) and < /i ^ 1 which finishes the proof. □ 
To conclude the study of r, we finally consider the case for large z. 



. 3 



Lemma 3.4. Let 5,oiq > be sufficiently small. Then the function f := 3|t|2, where r is defined 
in (j3.4|) . satisfies the bounds 

\di,d^,f{z,a^;h)\<Ck,ez'-'^+'' 
for all z £ {1 + |,a"^yo); a G (0,ao), < /i < 1 and k,i £ Nq. 



Proof. As in the proof of Lemma 13.31 we assume that 6,aQ > are so small that Lemma 13.21 holds 

3' ^ 3^ 



and zt{a^; /i) G (1 - 1, 1 + I) for ah a G (0, ao), < According to Eq. ([3^ . we have 



f{z,a^;h) = ^T{z,a'^;h)2 = / ^Ji- ^ e{a'^v?] h)du 

J zt{a'^;h) 



for z G (1 + 2' Q; ^yo)- Again, we split the integral as 



+ / =: A{a^;h) + B{z,a^;h) 



and, noting that 



we infer from Lemma 13.21 the bounds 



\d'^2A{a^;h)\ <C,\l + l-zt{a^;H)\^~' <Ce, £ G No 

since by the assumptions on oq and 6, we have |1 + | — zt{a'^;h)\ > 6 for all a G (0,ao) and 
< h <^ 1. In order to obtain suitable bounds on B{z, c?\ h), note that, by Lemma |2.H e satisfies 

since \u\ < a^^. Consequently, 2£{c^^ u'^ ] K) = {v?Ye^^\a'^v?;fi) behaves like a symbol and we 
infer the bounds 

\d'^2d':^,e{a'u'-M < Cfe,,(n2)^-^ fc,^ G No 
for M G (1 + |,a^-'^?/o)- This implies 

K^dti^ - ^ + ^(«'^'; ^)]l < Ck,iu-'^+^' 

and, since 1 — ^ + e(Q^?i^; /l) > 1 in the parameter regime we are considering, we obtain 

for all z G (1 + |,a~^yo)) ct G (0, ao); < ^ 1 and A;,£ G No, cf. Lemma fA.SI Trivially, the same 
bounds hold for A{a^] h) and the claim follows. □ 

3.2. Application of the Liouville Green transform. We have collected enough information 
on the function r to study the diffeomorphism ip. 



Proposition 3.5. Let C,, r he as in Lemma \3.1\ and Eq. (13. 4p . respectively, and suppose ao > Q is 
sufficiently small. Then the function Lp := C,~^ o r has the following properties: 

(1) For any a G (0, oq) and < ?i ^ 1, ip{-, a'^;h) : (0, a^^yo) — (0, Wo{a; h)) is an orientation- 
preserving diffeomorphism where wq is defined in (|3.3p . 

(2) ip(zt{a'^;h),a'^;h) = 1, i.e., ip maps the turning point to 1, 

(3) (fi is of the form 

(p{z, a^]h) = z[l + o?z^e2{z^ K)\ 
where the function £2 satisfies the hounds 

\di,d%{z,a^-h)\<Ck,fX^)-''^^' 

for all z G (0, a~^yo); « £ (0, ao); < ?i < 1 and k,£ e No, 

(4) ip{z, a^;h)^z and <p\z, a^; ^) ~ 1 for all z G (0, tt;o(a^; hi)), a G (0, ao) and < ^ <C 1. 

Proof. Throughout the proof we assume that ao is so small that zt{a'^; h) G (|, |) for all a G (0, ao) 
and < h<^l. 

(1) As already noted, the function ip{-,a'^;h) is an orientation-preserving bijection between 
the intervals (0, a~^yo) and {0,wo{a'^;h)). It is also clear from Eq. (j3.2p that (p{-,a'^;h) is 
smooth on (0,Z((a^;^)) U (zt(a^; ^), a~^yo) for any a G (0, ao); < h <^ 1 and smoothness 
around the turning point zt{a'^;h) follows from Lemmas 13.11 and 13.21 

(2) This is an immediate consequence of Eq. ()3.2p . 

(3) We distinguish different cases and start with z close to zero, say, z G (0, |). Since we 
stay away from the turning point, we may drop the fractional powers and, writing f := 
— ||r|2, ( := — ||C|2, we obtain ip = of. Now note that ( maps (0, |) to (— oo,C(|)) 
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diffeomorphically. Thus, C, := expo(^ : (0, |) — )• (0,e^''3)) difFeomorphically and, according 
to Lemma 13.11 C, is of the form 

^^^^ _ ^(ix) _ glogx-log2+l-x2£o(a;) _ Irf-^^-x^eoix) 

which imphes the bounds \C^''\x)\ < Ckx'^'^^^^^'''^^ for aU x e (0, |) and k £ Nq. Since 
(' ~ 1, the inverse satisfies the bounds \iC^^)^''\x)\ < CfcX™^^'f^~^'°^ for x £ (0,e^'^3)) 
and we note that = C^^ ° log. Now consider f := expof : (0, ^) — (0, '^)) which, 

according to Lemma 13.31 looks like 

and we obtain the bounds |9^25^f(z, q^; h)\ < Cfc,^z'"^^^i-'='°> for all z G (0, i), a G (0,ao), 
< /i < 1 and k,£ e Nq. Note that 

(f = C,~^ o f = o log o exp of = o f 

and this yields the bounds 

for all z G (0, ^), a G (0, oq), < ?i ^ 1 and k,£ Nq. As a consequence, since (p{z, 0; ^) = 
the fundamental theorem of calculus yields 

Lp{z , ; h) = z + / dpLp{z, I3;h)dj3 (3.5) 
Jo 

and we obtain ip[z,a'^;H) = z[l + a'^e{z,a'^;h)] where 

\di2d^,eiz,a^;h)\<Ck,i 

for all z G (0, |), a G (0, qq), < ^ 1 and k,£ G Nq. In order to prove that in fact 
e(z,a'^;h) = O(z^), we insert ip{z,a'^;h) = z[l + a'^e{z,a'^;h)] into Eq. (13. Ih which yields 

e'(0,a2;;i)[l + a^e(0,a^;n)] = 

and we infer that e'{0,a'^;h) = since e(0, q^;/i) = — ^ contradicts the above bounds 
on e. However, this and Eq. p.ip necessarily imply e(0,a^;/i) = 0. Consequently, the 
above bounds on e and the fundamental theorem of calculus show that in fact f{z, a^]fi) = 
z[l + z^£2{z,a^;fi)\ where £2 satisfies the same bounds as e. This settles the claim for 
(0,i). 

Around the turning point, i.e., for z G (|,3), say, it follows from Lemmas 13.11 13.21 13.31 
and [331 that 

\dl,d';^{z,a^-h)\<Ck^i 

for all a G (0, ao)) < ^ 1, /c, ^ G Nq and thus, the stated bounds are a consequence of 
the formula (13. 5p which is valid for all z G (0, a~^yo)- 

Finally, for the large z behavior, say, z G {2,a~^yo), we again remove the fractional 
powers and use the representation ip = o f where this time C, := ||C|^ and f := ||t|2. 
From Lemma l3.ll and Appendix |A] it follows that satisfies the bounds |(C~"'^)^'^H^)I — 
CkX^~^ for all X > Co > and /c G Nq where cq is arbitrary but fixed. Furthermore, from 
Lemma [3^ we infer that |5q2(/?(z, a^; < and thus, Eq. (13. 5j) implies 

ip{z,a^;H) = z[l + 0{z'^a^)]. 

The stated derivative bounds follow from the bounds in Lemma 13.41 and Eq. (13. 5p . 
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(4) For small the estimate <f'{z,a'^\ /i) ~ 1 follows immediately from the above and for large 
2 it is a consequence of 93 = C,~^of, Lemma [3. 11 and the fact that f'{z, a^; h) ~ 1, uniformly 
in small a, h, which can be read off from the definition of r. The estimate ip{z,a^;h) ~ z 
now follows from ip{z, a^; h) = Jq (p'{u, a^; h)du. 

□ 

Remark 3.6. From Proposition 13.51 we can also read off the behavior of WQ{a'^;fi) as a — )• 0+. 
Indeed, we have 

wo{a'^;h)= lim (p{z,a^\h) lim z ^ a^^ , 

uniformly in small h. 

Now we can apply the Liouville-Green transform induced by the diffeomorphism ip{-,a'^;h) to 
Eq. (USD. 

Proposition 3.7 (Normal form reduction). Assume 0<a<^l,0<h<^l and let ip(-,a'^;h) be 

the diffeomorphism from Proposition \3.5[ Furthermore, as before, set = ^ + 4£'^ . Then the 
function g, defined by 

9 (<^(ie-t,a2.;,)^ = e-t(/?'(fe-t,a2;;i)f/(x), 

satisfies the equation 

- hlg"{w) + (1 - ^) g{w) = h^Vi{w, K)g{w), hi = ^ (3.6) 

for w G (0, tt;o(a^; /i)) if and only if f is a solution to Eq. (j2.3p for x > xq. Here, V\ is a suitable 
function that satisfies the bounds 

\di2dtViiw,a^;h)\ < CkiH-''+^', k,£ G No 

in the above domain ofw,a,h. Furthermore, there exists a unique solution g^(-,E;h) to Eq. (j3.6p 
with the asymptotic behavior 

g+{w,E;h)r^(^^y~''^ {w^0+) 

and the outgoing Jost function /+(•, E; h) of the operator H{h) is given by 

f+{x,E-h) = e^(p'{^e-^,a^;h)-^g+ (^^{^e'^ ,a^;h), E;h^ . 

Remark 3.8. The crucial point of Proposition 13. 71 is that the right-hand side of Eq. (j3.6p is globally 
small, even after dividing by h-f. 

Proof of Proposition \3. 7| By construction of ip{-,a'^;h) (see Eqs. (|3.ip and (|2.2p ). the function 
g{ip{z,a'^;h)) = ip'{z,a'^;h)2g(z) solves Eq. (|3.6p on w ^ {0,wo{a'^;h)) with 

2..^ 2.,^_ ^? ^3(^"(z,a2;n)2 l^"'{z,a'-h) 



Vi{^iz,a'; n),a'; h) = i - J \' - ^ ' , (3.7) 

\4 (/3'(z, a^; n)^ 2 ip'{z,a'^;n)'^ J 

if and only if g solves Eq. (j2.6p on 2: G (0, yo) where yo = 2e~ 2 . However, this is the case if and 
only if f{y) = g{^) satisfies Eq. (j2.5p on y G (0, yo) which, according to Lemma [2Tl is true if and 
only if f{x) = e^f{2e~^) solves Eq. (j2.3p on x > xq. The asymptotic behavior of 5+ follows from 
ip{z,a'^]h) ~ z and (p'{z,a'^;h) ~ 1 as z — )■ 0+, see Proposition 13.51 which implies the asymptotics 
ip~^{w, o?] K) ^ w and ((^~-^)'(r«, ; ^) ~ 1 as — )• 0+ for the inverse V3~"^(-, o?\ h) of (/?(•, a^; H). 

It remains to prove the bounds on Vi. To this end, note first that, according to Proposition 13.51 
we have 
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and, similarly, 



on z G (0, a^-'^yo) since Lp'{z^o?\K) ~ 1 and wi^{q?\K) ~ a^^ ^ see Remark 13.61 This proves 
\y\{w^o?\fi)\ < 1 in the present parameter regime. Now note that Proposition 13.51 implies the 
bounds 

for all z G (0, a^^yo) and we infer inductively from 
the bounds 

\di,dt^-\w,a^;h)\ < CkAw)'''^^\ k,£ G No 
for all w G {0,wo{a'^; h)), cf. Lemma [A. 51 below. Consequently, by setting z = ip~^{w,a^;h) in 
Eq. (|3.7p . the claimed bounds on Vi follow via the chain rule and the bounds on ip from Proposi- 
tion [331 cf. Lemma |A.3[ □ 

4. Construction of a fundamental system for the normal form equation 



In this section we construct a fundamental system to Eq. (j3.6p . This is done by perturbing 
modified Bessel functions which are solutions to the "homogeneous" version (i.e., right-hand side 
set to zero) of Eq. (|3.6p . The fundamental system is constructed by solving suitable Volterra 
equations the kernels of which are composed of modified Bessel functions. This requires a good 
understanding of Bessel functions and we refer the reader to Appendix [D] for the necessary results. 
The main technical difficulty we are faced with comes from the fact that we need good estimates for 
the errors and all their derivatives. In order to keep things as readable as possible, we distinguish 
different regimes, namely exponential / oscillatory as well as j- small / large. In fact, we proceed 
analogously to the construction of the Bessel functions in Appendix [Dl 

4.1. Preliminaries. We start with an elementary result that proves to be very useful in the sequel. 

Lemma 4.1. Let I cM. be an open interval and consider the inhomogeneous differential equation 

u"{x) + V{x)u{x) = fix) (4.1) 

for X € I with u, V, f complex-valued and V, f continuous on I. Suppose further there exists a 
function uq that does not vanish on I and solves the homogeneous equation 

u'^{x) + V{x)uo{x) = 0. 

Define 



a{x) :- 



Uo{x')f{x')dx' 



Uq (X )ax 

where xq £ I is an arbitrary constant. Then u{x) := Uq{x)[1 + a{x)] is a solution of Eq. (|4.ip . 

Proof. This can be verified by straightforward differentiation. □ 

In the following, the ratio 2-^ will appear frequently as a parameter and therefore we use the 
abbreviation u := 2^. The various parameters and their dependencies are summarized in Table [H 



4.2. A fundamental system of the normal form equation for small i/. We start the con- 
struction with the case < < 1 which means E <h (and also hi ~ 1). Thus, the classical Bessel 
asymptotics Eqs. (1D.2P and (ID.4P become relevant. The first result deals with the exponential 
regime, i.e., to the right of the turning point and towards infinity. 
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Table 1. Parameters 



Parameter Relevant range Useful relations 



E 


< S < 1 




h 


< ;k 1 






< a < 1 






< ^1 < 2 




-=2f 


< < oo 





4.2.1. T/ie exponential regime. After rescaling by hi = — , i.e., setting G{v) := g{hiv), the "homo- 



geneous" version of Eq. (j3.6p takes the form 

G"{v) + ( - 1 ) G(^) = (4.2) 



1.2 + 1 



for u = 2^ which is a modified Bessel equation. Eq. ()4.2p has a fundamental system {Bj(-, u) : j = 
1,2} of the form 

Bi{v,iy) = e--[l + bi{v,iy)] 
B2{v,u) = e^[l + b2{v,u)] 

where the functions bj{-,iy) are real-valued and satisfy the estimates 

\did^bjiv,,^)\ < C,,i{v)-'-\ kjGNo 

for all t> > 1, say, and < < z^o (where z^o > is some fixed constant). This follows by appropriate 
Volterra iterations and is explicitly proved in Appendix [Dl see in particular Lemma lD.4[ 

Lemma 4.2. Let < E < h 1 and wi > be sufficiently large. Then there exists a fundamental 
system {gj{-,E; h) : j = 1, 2} for Eq. (|3.6|) on [wi,WQ{a'^; h)] of the form 

gj{w, E; h)=Bj{^, iy)[l + hcj{w, E; h)] 
where a, v are given in TableUi The error terms Cj{-,E;h) are real-valued and satisfy the bounds 

\d'Edtcj{w,E;h)\ < Ck,eh-', k,ieNo 
for all w G {wi,WQ{a^;h)\, and in the above range of E, h. 

Proof. We start with the growing solution g2{-,E;h) and according to Lemma l4. II we consider the 
equation 

hc2{w,E;h) = - / B2{^,iy)-'^duB2i^,iy)'^a'^Vi{v,a'^;h)[l + hc2iv,E-h)]dv 

Jwi Jv 

for w G [wi,WQ{a'^; h)] which is well-defined since B2{^,iy) > for all w > wi provided that wi > 
is sufficiently large (recall that /ii ~ 1 in the domain of E and h which is being considered here). 
Rescaling by hi yields 

/•X i-x 

hc2ix,E;h) = - / B2iu,i^)-^duB2iy,i^fh'^Viihiy,a^;h)[l + hc2iy,E;h)]dy 

for X £ [h'^^wi, h^^WQ^a'^; h)] where 62(10, E; h) = C2{^,E; h). This equation is of the form 

fX j-X 

hc2{x,E;h)= / e''^''a{u,E;h)due^yb{y,E;h)[l + hc2{y,E]h)]dy 
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for suitable functions a, b satisfying 

\d'Ed^Mx,E;h)\ < CkA^)-^^-', \d'Ed^,bix,E;h)\ < CkAxr'^h''-' 

for all X G [h^^wi,h^^wo{a'^;h)], < E < h 1 and k,£ E Nq. Furthermore, we recall the 
bounds \d^h^^wi\ < Cfh~^ and |(9^/ij~"^it;o(a^; K)\ < C^h^'^^^ in the relevant domain. Consequently, 
Proposition IB. 21 yields the claim concerning g2{-,E;h). 

The second solution gi{-,E]h) can now be obtained by the standard reduction ansatz. More 
precisely, we set 

/"Wo{a'^;h) 

gi{w,E;h) = -2g2{w,E;h) / g2{v, E;hy'^dv 

J w 

for w € [wi,wo{a'^;h)]. This is certainly well-defined provided that twi > is sufficiently large and 
h > is sufficiently small. It is straightforward to check that gi{-,E; h) is indeed of the stated form 
(cf. Lemma EI]). □ 

4.2.2. The oscillatory regime. Next, we construct a fundamental system in the oscillatory regime. 
To this end note that it follows by Frobenius' method that Eq. (j4.2p has solutions B±[-,u) of the 
form 

B±{v,u) = v^^'''[l + h±{v,u)\ 

where the error terms satisfy \dld^b±{v, u)\ < Ck,iv'^^^^'^~^'^^ for aU v G [0, vq], u ^ (0, z^o]; A;, ^ G No 
and some fixed vq,vq > 0, cf. Corollary ID. 21 Furthermore, it is known that \B±{v,i')\ does not 
vanish for any v > 0, see Corollary ID. 21 Consequently, we construct a fundamental system of 
Eq. dM]) by perturbing B±{^,u). 

Lemma 4.3. Let < E < h <^ 1 and wi > 0. Then there exists a fundamental system {g±{-, E; h)} 
for Eq. (j3.6p on (0, wi] which is of the form 

g±{w,E;h)=B±{!^, y) [1 + h^c± {w, E; H)] 
where, as before, = ^ + AE"^ and v = 2-|. The error terms c±{-,E; h) satisfy the bounds 

\d'Edlc±{w,E-h)\ < Ck,ew^''h-', k,£eno 
for w G (0,tt;i] and in the above range of E, h. Finally, we have c±{0,E;h) = dwC±{0,E;h) = 0. 
Proof. According to Lemma |4.1| we have to construct a solution to 

h^c±{w,E;h) = - / B±{^,u)-^du B±{f,iyfa'^Vi{v,a'^;h)[l + h'^c±{v,E;h)]dv (4.3) 

Jo Jv ^ ^ 

with hi = ^. By introducing the new variable w = hie~^ we obtain 

poo py 

h^c±{x,E;h) = - / e-^B±{e-\u)-^dze-yB±{e-y,vf 

J X J X 

X h^Vi{hie-y,a^]h)[l + h^c±{y,E;h)]dy 

for X G [— log^,oo) where c±{w,E;h) = c±{— log E;h). Now recall that B±{e~^,i') = 
e~^2='=*^)^[l + b±{e~^ , v)] and trivially, we have the bounds 

\d'Ed',e-^[l + b±{e-\u)]\<CkX^r''f^-', kjeno 

by noting that (?£;z/ = 2h~^. Furthermore, according to Proposition 13.71 we have the bounds 

\d'Ed^h'^e-''Vi{hie-'',a^;h)\ < Ck/{x)-'^-^h^-\ A;,^ G No 
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in the relevant domain by noting that a > h. This shows that the Volterra equation for c± is of 
the form 

J X J X 

where the functions a, b satisfy the estimates 

for ah E Nq and in the relevant domain of x, E, h. As a consequence, we obtain from 
Proposition IB.ll the existence of c± with the bounds 

\d'Ed',c±{x,E;h)\ < C,Ax)-'-'h'' 

for all x £ [— log oo), < E < h <^ 1, k,i G Nq and via the chain rule, these bounds translate 
into the claimed ones for c±. Finally, the fact that c±(0, E; h) = 5^c-i-(0, E;h) =0 follows directly 
from Eq. □ 

Remark 4.4. The bounds on c± in Lemma 14.31 are not optimal. In fact, it is not hard to see that 
one has 

\d'Edic±{w,E-h)\ < Ck/w'^''''^'^'^'^^h-^ 
but we leave it to the interested reader to prove these stronger bounds. 

4.2.3. Matching of the fundamental systems. As a next step, we glue together the two funda- 
mental systems obtained in Lemmas 14.31 and 14.21 This amounts to calculating the Wronskians 
i?; (•, -E; /l)) or the asymptotics of g±{w, E; h) for large w. 

Lemma 4.5. Let wi > be sufficiently large. Then the functions g±{-, E; h) from Lemma \4 . 3| have 
the representation 

g±{w,E;h) =TT-2 2-2^"'T{l±iu) e—[l + h^2{E;h)][l + hd2{w,E;h)] 



aw , 



- ie^^^e" — [1 + h^i{E; h)] [1 + hdi {w, E; h)] 

where the error terms dj{-,E;h), j = 1,2, are real-valued and we have the bounds 

\d'Edtdj{w, E-h)\< Ck,ih~', Id'^jjiE; h)\ < c,h-' 
for all w G [wi,wo{a'^;h)], < E < h <^ 1 and k,£ e Nq. 

Proof. Since {g±{-,E; h)} and {gj{-,E] h) : j = 1, 2} are fundamental systems for the same equation, 
there exist connection coefficients "y±j{E; h) such that 

9±{w,E;H) = -f±A^;H)gi{w,E;h) + -/±A^; h)g2{w, E; h) 

and this implies 

^._ W{g±{;E;h),g2{;E;h)) W{g±{; E;h), g^j-, E;h)) 

^^•'^ ' ' W{gi{;E-h),g2i;E;h))' ^^''^ ' > W{gi{; E; h), g2{; E;h)) ' 
We have 

Wigi{;E; n),g2i;E; h)) = lW{B^{-, z.), i?2(-, + 0{h)] = |-[1 + 0{h)] 

where the 0~term has the property that d^EO{h) = 0{h^^^) in the domain under consideration. 
Analogously, we have 

W{g±{; E; h),g,{; E; h)) = lW{B±{-, u),B,{; i^W + Ocih)] 

and in particular, 

W{g±{; E- h),gi (•, E; h)) = -/i^ V-525±-r(l ± iv)[l + Oc{h)l 
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see Lemma lD.3[ This shows 7±^2(-E'; /i) = vr 22 2 ='=*'^r(l it zi/)[l + Oc(^)] and the claim follows. □ 

Based on Lemmas 14.21 14.31 and 14.51 we obtain the desired representation of the semiclassical Jost 
function of the operator H{h). 

Corollary 4.6. Let < E < h <^ 1. Then the semiclassical Jost function /+(•, E; h) of the operator 
H{K) has the representation 

E- h) = 2-i+-/i^-ef (/.'(|e-t , a^; h)-hg_{^{le-^^ ,a^-h), E; h) 

for all X > xq where, as always, = + iE"^ , v = 2-^ and the function g-{-,E; h) is given in 
Lemma \4.3\ 

Proof. According to Proposition 13.71 the solution g^{-,E; h) of Eq. (j3.6p . which corresponds to the 
semiclassical Jost function, has the asymptotics 

g+{w,E;h) ^2-^+'''{aw)'2-''' {w^0+). 
On the other hand, the solutions g±{w,E; h) from Lemma 14.31 satisfy 

g±{w, E; h) ~ B±{^^) ~ h^'^^'" {awf^^'" {w ^ 0+), 

see Corollary Eil and this implies g+{-,E;h) = 2-^+'''h^-'''g_{-,E;h). □ 

4.3. A fundamental system of the normal form equation for large u. In this section we 
study the large case, i.e., E ^ h. It is important to note that in this regime the parameter 
hi = ^ = {y^ + 3)^^ can become arbitrarily small. The construction relies on the large parameter 
asymptotics of the Bessel functions based on Airy functions, see Appendix |Dj Consequently, it is 
convenient to transform Eq. (j3.6p according to Appendix [Dj First, we rescale and introduce the 
new independent variable v := Setting G{v) = g{hiuv), Eq. (|3.6p transforms into 

G"{v) - u\l - ^)G{v) + ^G{v) = -a^Viiv, E; h)G{v) (4.4) 

where v G (0, {hiv)^'^WQ{a^;h)) and Vi{v,E;h) := {hiv)'^Vi{hivv,a^\h). At this point is useful to 
note that hiv ~ 1 if > 1 and thus, we have rescaled by a harmless factor. From Proposition 13.71 
we have the bounds 

\d'Ed';Vi{v,E-h)\<CkA^)-''h~' 

in the respective domain of f , E and h. Following Appendix [D] we now apply the Liouville-Green 
transform 

Hav)) = C'{vf'G{v) 
with C, : (0, 00) — )• M from Lemma 13.11 From Eq. (j4.4p we obtain 

^"{0 - ^\H0 - ^2(C)<f (C) = -«V3(C, E- hMC) (4.5) 

for C G (— 00, Co(-S; /i)) where Co{E;h) := C{{hii^)~^wo{a'^;h)) ~ and 

VsiCiv), E;h) :=C'{vr^Vi{v,E;h)=av)^^Vi{v,E;h). (4.6) 



Here we use the symbol to denote both, the diffeomorphism from Lemma 13.11 as weU as the new independent 
variable. 
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The function V2 is given explicitly in Appendix [D] and satisfies jVg (C)| ^ C^iC) ^ ^ for all C G M 
and k gNq. Note further that we have the bounds 

\d'Ed'^Vs{C,E;h)\ < CkACr^-'h-', c < o 
Id^jd'^Vsic, E;h)\ < CkAO'-'h-', c > o 

from Proposition 13.71 and Eq. (|4.6p . 

4.3.1. The exponential regime. It is a standard fact that the homogeneous version of Eq. (|4.5p is 
solved by Bessel functions, see Appendix [Dj In particular, for C ^ 0, there exists a fundamental 
system v) : j = 1, 2} of 

$"(C) - y\HC) - V2{CmC) = (4.7) 

which is of the form 

</.i(C,z^) = Ai(z/ic)[i + 1^-^01 (C,/^)] 

</.2(C,z^) =Bi(z.ic)[l + i^-^a2(C,z^)] 

where \aj{C,,v)\ < 1 and 

\did'^aj{C,u)\<Ck^e{0-l-''iy-', C>1 

as well as 

|a5a,(0,i/)| < C,l/-^ \did^aj{0,u)\ < dJ-' 
for all 3> 1 and k,i £ Nq, see Lemma [D.4[ Based on the Bessel functions (f>j{-, i') we now construct 
a fundamental system of Eq. (j4.5p in the exponential regime. 

Lemma 4.7. For G [0,Co(^i^)] there exists a fundamental system {^i{-,E;h),^2{'TE;h)} for 
Eq. (j4.5p of the form 

^i{C,E;h) = MC,^)[^ + fiM(,E;h)] 

^2{C,E;h) = M(,^)[^ + h<^2{C,E;h)] 

with V = 2-^. The functions aj, j = 1, 2, are real-valued and \o'j{C, E;h)\ < 1 for G [0, Co{E; h)], 
< h <^ E <^ 1. Furthermore, they satisfy the bounds 

\d'Ed'^a,iC,E;h)\ < Ck^eicy^h'', C e [hCo{E;h)] 

as well as 

|a^aj(0, E-h)\< Cih-\ \d^Ed^aj{0, E;h)\ < Civlh'^ 
for alio < h<^ E and k,£ £ Nq. 

Proof. We start with the growing solution and consider the Volterra equation 

ha2{C, E; h) = -0? <P2{r]\ ^)'^dr]' ^2iv, 1^)^3(7?, E; h)[l + ^1^2(77, E; h)]dr] 

Jo Jt) 

2 

for C G [OXo{E;h)]. Rescaling by z^~3 yields 

a2{x,E;h) = -a^u-^ / / (l)2{v'^^u,v)-'^du(t)2{v-iy,ufVi{v~^y,E;h)[l+a2{y,E;h)]dy 

Jo Jy 

for X G [0,u^Co{E;h)] where a2{x,E]h) = ha2{v~^x,E;h). We have the bounds \d^EV^C,Q{E;h)\ < 
Cihrs~^. From Eq. ()4.6p and Proposition 13.71 we infer the estimates 

\d'Edyiy-"sV^{u-ly,E-h)\ < CkAv)'^^^-' 

^In fact, V3 decays much faster as ^ — > —00, namely exponentially. However, we do not exploit this fact. 
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in the relevant domain of y,E,h and for all k,i £ Nq. Thus, we obtain from Proposition IC.8I the 
existence of a2{-,E;h) with 

\dj,d^,a2{x,E;h)\ < CkAx)"2-^m-' < Ck,,{x)-''h^-^ 

and the stated bounds for ai follow. 

The solution E; K) can now be constructed by the standard reduction ansatz, i.e., by setting 

The so-defined ^i{-,E; h) is clearly a decaying solution to Eq. (|4.5|) and by using VF(Ai, Bi) = 
see e.g., [H], it is not hard to see that ^i{-,E;h) is indeed of the stated form. □ 

4.3.2. The oscillatory regime. Similarly, in the oscillatory regime C ^ there exists a fundamental 
system {4>±{-,i^)} of Bessel functions for Eq. (|4.7p of the form 

MC, = [Ai(i^ic) ± iBi(i.ic)][l + y~^a±{C, y)] 

where the error terms a± satisfy analogous bounds as Oj above in the respective domain of C, and v, 
see Lemma lD.51 By perturbing the basis {</'±(-, v)} we obtain a fundamental system for Eq. (j4.5p . 

Lemma 4.8. For C ^ there exists a fundamental system {(^±[-^E] /i)} of Eq. ()4.5p of the form 

^±{CE-h)=(l)±{C,y)[l + ha±{(:,E-h)] 

where v = 2-^ and \o'±{(, E;h)\ < 1 for C ^ and < h E 1. Furthermore, the error terms 
a± satisfy the bounds 

as well as 

\d^Ea±{0,E;h)\ < Ceh'^ \dljd^a±{0,E;n)\ < Cei^h'^ 
for all < h <^ E <^ 1 and k,£ £ Nq. 

Proof. By choosing sufficiently large, we obtain from Lemma [D. 5 1 that \<j)±{C, z^)| > for all C ^ 0. 
Thus, the function cr_ can be constructed by solving the Volterra equation 

ha^{C,E-h) = -a^ (t>^{'n',uy^dr]' <P^{ri,ufV3{i],E;h)[l+ha^{ri,E;h)]dri. 

J — oo J q 

This is done via Propositions IC.7I and IC.81 completely analogous to Lemma 14.71 (or Lemma ID.SP . 
is just the complex conjugate of □ 

4.3.3. Matching of the fundamental systems. Next, we glue together the two fundamental systems 
{<I>j(-, E] h) : j = 1, 2} and E; h)}. To this end we make use of the representation 

for C > 0, 1 where \dl,a±j{v)\ < Civ~^ for all £ G No and j = 1, 2. This representation follows 
easily from the fact that the Airy functions are defined globally, see Lemma ID. 61 for an explicit 
proof. 

Lemma 4.9. For C > the functions ^±{-,E; h) from Lemma \4.8\ have the representation 

^±{C, E; h) =[l + v-^p±,i{E; n)]cI>i(C, E- h) ± i[l + u'^f3±,2{E; h)]<^2{C, E; H) 

where E;h), j = 1,2, are from Lemma \4. 7| and v = 2j-. The functions f3±j satisfy the 
estimates 

\d'E(3±AE;h)\<c,h-' 

for all < h -t: E 1 and £ € Nq. 
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Proof. By construction (see Lemma \T7l^ we have W{^i{-, E; h),^2{-, E; h)) = z^s VF(Ai, Bi). Fur- 
thermore, by evaluating the Wronskians at = and noting that h < z^^^, we obtain from Lemmas 
imOl that 

W{^±{;E;h),<^j{;E;h)) = Jw{A±,Aj)[l + u-^P±j{E;h)], i = l,2 

for functions j3±j{E\K) satisfying \d^^j3±j{E\h)\ < Cih~^, I G No, where A\ := Ai, A2 := Bi and 
A± := Ai lb iBi. The result now follows by a straightforward computation, cf. Lemma lD.61 □ 

Now we are ready to obtain the desired representation of the outgoing Jost function in the case 
of large ly. 

Lemma 4.10. For < /i <C -E <C 1 the semiclassical outgoing Jost function f^{-,E;h) of the 
operator H{H) has the representation 

f+{x,E;h) =vr52-ie^(^+f)?i5--i.i— et(^'(|e-§,a2;?i)-K'(i;^(ie-i,a2;?,))-| 

x^^{a^^y,{^e-^,a';h)),E;h) 

where v = 2j-, = ^ + 4£'^ and <!*-(•, E] h) is given in Lemma \4-8[ 

Proof. By construction, the functions g±{-,E; h) defined by 

g±{w,E-h) ■.= C'{^)-U±{Q{^^lE-h) 

with $-!-(•,£'; h) from Lemma [4.81 form a fundamental system for the normal form equation (|3.6p on 
(0, 1] provided that v = 2j- 1. According to Lemmas l4.8llD.5l and the standard Airy asymptotics 
(see Corollary IC.4P we have 

^±{C„E;h) ~ 0±(C,i^) ~ Ai(i/ic) ±iBi(z^ic) ~ 7r~^e=^*f (-i/ic)"^e=F^i''(-^)''" 
as C — ^ — 00 and thus, 

g±{w, E; h) ~ vr-ie±*f z.-ic'(^)-i [-C{^)]-\e^"'"^^-^^^^^^"" 
as w ^ 0+. Now recall from Lemma |3. II that 

I [-C(x)] i = - log X + log 2 - 1 + 0(x2) 
C'(x)[-C(x)]5 = i + 0(x) 

— )• 0+ which imply 

~ / 7^ t^ -- ±i- --I w \ \ T«;^[-log(Tr^)+log2-l] 
g±{w, E;h) TT 2e^^4,l/ 6(^j2g^ i s^^^^; S J 

r^7r-h'^'^e^'^''+^^h-^^'''iy~ri^'''{aw)^^''' 

as U) — )■ 0+. According to Proposition 13. 7|, the solution g^{-,E;h) of Eq. (|3.6|) which corresponds 
to the outgoing Jost function /+(•, E; h) of the operator H{h) has the asymptotics 

g+{w,E;h) r^2'^+'''{aw)'2~^'' {w ^ 0+) 

and thus, we obtain 

g+{w,E;h) = TTh-h''^''+^'>h^-'''J-'''g^{w,E;h) 

which yields the claim by undoing the transformations that led to Eq. ()3.6p . □ 
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5. The scattering matrix and the spectral measure 

In this section we compute the scattering matrix as well as the semiclassical spectral measure 
associated to H{h). Thereby, we complete the proofs of Theorems ll.2l and ll.3i The most important 
step consists of computing the Wronskian W{f-{-, E; h), f^{-,E] h)). To this end we need a similar 
asymptotic description of f-{-,E]h) as we have obtained in the previous section for f+{-,E]h). 
This, however, is trivial due to the assumed symmetry of the problem. In fact, the corresponding 
result for f-{-,E;h) can be obtained from f^{-,E;fi) by switching from x to —x. From now 
on we write 99+ instead of 99 for the diffeomorphism in Proposition 13.51 and ip- will denote the 
corresponding diffeomorphism associated with the construction of f^{-,E;h). Furthermore, it is 
useful to introduce the following new notation. 

Definition 5.1. By e{x,E;h) we denote a generic real-valued function that satisfies the bounds 

\d'Ed^,e{x,E-h)\<Ck,ih-''-' 

for all k,£ G Nq and in a domain of x,E,h that follows from the context. Moreover, we write 
ec{x,E; h) if the function attains complex values as well. 

5.1. The Wronskian. Based on the results of Sec. H] we have a good description of f^{-,E;h) 
on (xo,oo) and similarly for f-{-,E;h) on (— 00,— xq). Without loss of generality we assume that 
xq < and compute the Wronskian at x = 0. 

Lemma 5.2. For the outgoing semiclassical Jost function /+(0, E; h) of the operator H{h) we have 
the representation 

/+(0, E; h) = a^-f+iE; fi)eTi'^s+(E;n)+iT+{E;h)) N ^ f^^^-^. ^ ^y-ls+{E-fi) 
/;(0, E- h) = ah~^c+{E; h)-/+{E; n)e\iS+{E-fi)+iT+{E-fi)) U ^ j^^^^. ^ ^y-ls+(E;h) 

where C-^{E;h), S^{E;h) > 1, T^{E;h) S M, \ j^{E;h)\ ~ 1, and we have the bounds 

\dj,A{E;h)\<Ceh-', A G {c+,^+, S+,T+} 
for all £ eNo, < E -^1 andOKh-^l. 

Proof. We start with the case < < 1, i.e., a ^ h. According to Corollarv 14.61 and Proposition 
we have 

f+{x,E;h) =i7r-U5-^^r(l -iz.)[l + n72(^;^)]etv?V(|e-t,a2;;i)-|efv+(|^"*."';^) 
"1 + he{x, E; h) + ee(x, E; ;i)e-2t^+(f , 

cf. Definition 15. H where, as always, = ^ + AE"^ and v = 2j-. In order to compute f^{0,E;h) 
we note that 

V ' 

0(a) 

= e{E-h) 
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by Proposition 13.51 Hence, we infer 

/;(0,ii;;;i) = i7r-U-^-r(l-iz.)[l+72(i?;;i)]v^',(^,a2;^)lef'^+(^'°'^^) 

By setting S-^{E;h) := aip+{^,a'^;h) we obtain from Proposition 13.51 the bounds \d^^S+{E;h)\ < 
Ceh~\ £ G No as well as S+{E;h) > 1. Furthermore, we set T+{E;h) := -2£;iog^ and the 
result follows by recalling that ip'^{^, a^;h) ~ 1 and /i ~ a in the parameter regime which we are 
considering. 

Next, we turn to the case v ^ 1, i.e., E ^ h. From Lemmas 14.101 and 14.91 we have the 
representation 

f+{x,E;h) =^^2-5e^(-+f)n^-z.§— ef^V(|e-f,a2;;,)-iC'(^(^+(^e-§,a2;;,))-^ 
"[1 + i^-'P-,i{E; h)]^^iC{^^+{^e-^,a^; h)),E; h) 

- i[l + v~^P-,2{E; K)]^2{at^+{h-^^.c^^; ^)), E; h) . 
Now recall from Lemma 14.71 (see also Lemma ID. 41 and Corollary IC.2j) that 
^i{C,E-h) = <l,i{C,u)[l + hai[Q,E-h)] 

= Ai(i/§C)[1 + i^"'ai(C, + ^^i(C, E- h)] 
= (47r)-|(i.ic>-^e-4C^/^[l + z.~iai(C, + haiiCE; h)] 
for C G [OXo{E;h)] and analogously, 

^2{C, E; h) = (z.fo-ie'^i^''' [1 + v-\2{C, ^)] [1 + ?i^2(C, E; h)] 
for real-valued aj{-, E; h),aj{-, E; h) that satisfy 



for all C £ [1) Co{E; h)], k,£ £ Nq, j = 1,2 and in the relevant domain of E, h. Now note that for 
|x| ^ 1 we have 



and thus, 

\u-'a,{C{i;^+{^e-^^,a';h)),u)\<h. 
Consequently, f^{x,E;h) can be written as 

U{x, E;h)=- 2-^ie^(-+f et^V(le-f , a^; e"! , a^; 

x[l + u-'p^,2{E-,h)]{.H{^^^+{h-^,a';h)))-h'^l^^^^^^ 
1 + he{x, E; h) + edx, E; h)e-^''l'^^^'^+'^h'^ 



for I a; I < 1. By definition of C (see Lemma |3.1|) we have 

l^(x)i = x[i + o{x~^)], c'ixyhixy^ = 1 + o{x~^) 

for, say, x >2 and the O-terms behave like symbols. We therefore obtain 

C'(i;^+(|,a^;n))-|(.lc(i;^+(|,a^;n)))-ie4^(^--(i-^^^'))^^^ 
= u-^l + he{E;h)]e^^+^^'^'^ 
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where 



and by noting that hv c:^ a this yields the claimed expression for f^(0, E;h) with T^{E;h) = 
—2E log a. The expression for /+(0, E; H) follows immediately by observing that differentiation with 
respect to x only affects real-valued terms and by the symbol behavior of the involved quantities 
one picks up a factor h~^. □ 

Due to the symmetry of the problem we find the analogous representation of the Jost function 
f-{0, E;h) by simply replacing "+" in Lemma 15.21 by "— " and attaching a minus sign in front of 
the expression for the derivative. Consequently, we can easily calculate the Wronskian. 



Corollary 5.3. The Wronskians W{f-{-,E;h), U{-,E;h)) andW{f-{-,E;h), f+{-,E;h)) have the 
form 

W{f_{; E; h),U{; E; h)) = c{E; h)^_ {E; h)^+{E- ?i)an~iei(^(^'^')+^^(^''"')) [1 + he,{E- h)] 

Wif.{.,E;h),Ui.,E;h)) = ciE;h)^iIiE^h^iE;h)an~'e'^('(^^^^^^^^^ 

X [1 + hec{E; h)] 

where 5 = 5_ + 5+ , T = T^ + T^, c = c- + c^ and c±, 7± , S± as well as T± are from Lemma 15.^1 
Proof. From Lemma 15.21 and the above comments we obtain 

Wif.{;E; h),f+i;E; h)) =[c-{E; h) + c+{E- h)]j.{E; h)j+{E; h) 

X a;i-iei(^(^''"')+^^(^''"^))[l + he,{E- h)] 



for S = S- + S+,T = T^+T+ and similarly for W{f^{-, E; h),f+{-, E; h)). 



□ 



5.2. Proof of Theorem 11.21 Recall that the reflection and transmission amplitudes r(E;h) and 
t{E;h) are defined by the relation 



E; h) = t{E; h)f+{x, E; h) + r{E; h)f_{x, E; h). 

Consequently, we obtain 



W{f.{-, E- h),f_{-,E; h)) = t{E- h)W{f.{-,E; h), f+{-,E; h)) 



and the claimed expression for t{E;h) follows from Corollarv 15.31 Similarly, 



W{f.i-,E; h),f+i-, E; h)) = r{E; h)W{f.i-, E; h),f+{-, E; K)) 
and Corollarv 15.31 vields the claimed form oi r{E;h). 

5.3. Proof of Theorem 11.31 In order to prove Theorem 11.31 note that the semiclassical spectral 
measure e(0, 0, E; h) is given by 



e(0,0,£;;;i) = Im 



/_(0,£;;n)/+(0,i^;/i) 



Im 



W{f.{-,E;h)J+{;E;h)) 
According to Lemma 15.21 we have 
f\{Q,E-h) f'_{0,E;h) 
f+{0,E;h) f-{0,E;h) 

= [c+{E; h) + c-{E; h)]h-^ 1 + he{E; H) + ec{E; 7i)e" 
and hence, 



r+{0,E-h) f_iO,E;h) 
f+{0,E;h) f-{0,E;h) 



-1 



iS+{E;K) 



+ ec{E;h)e--^^-^^'^^ 



Im 



f.{0,E;h)f+{0,E;h) 
W{f.{-,E-h)J+{;E-h)) 



h 



c+{E;h) +C-{E;h) 
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e{E]h)e' 



\S+{E;h) 



which is the claim. 



Appendix A. Symbol calculus 



Roughly speaking, a function is said to have symbol character if it behaves like a polynomial 
under differentiation, i.e., each derivative loses one power. More precisely, we give the following 
definition. 

Definition A.l. A function / : I — )• M is said to belong to the set Sa{I) for a E M and / C M 
open iff / is smooth on / and satisfies 



for all A; G No and all a; G / where Ck > is a constant that only depends on k. Elements of Sa{I) 
are said to be of symbol type (or have symbol character, behave like symbols). 

The point is that symbol behavior is preserved under the usual algebraic and differential opera- 
tions. 

Lemma A. 2. If f e Sail) and g G Sa{I) then f + g e S^il), A/ G Sa{I) for any A G M, and 
f G Sa-i{I)- Furthermore, if f £ Sa{I) and g G 5^(1) then fg G 

Proof. The statement concerning the product fg is a consequence of the Leibniz rule and the rest 
follows directly from the definition. □ 

As a very convenient fact, the symbol behavior is even preserved under composition of functions. 

Lemma A. 3. Let f G Sail) and g G Sjs{J) with g{x) G / and \g{x)\ > \x\^ for all x £ J. Then 
f o g £ Sap{J). 

Proof. Since g has range in /, the composition f og \s well-defined and smooth on J. Furthermore, 
we have \ f{g{x))\ < [(^(x)!" < \x\°^^ for all x G J (use either \g{x)\ < \x\l^ or \g{x)\ > \x\^ depending 
on the sign of a). In order to estimate (/ o g)^''^ for A; G N, we claim that 



for suitable functions gj^k £ Sjp-kiJ)-, i = 1, 2, . . . , A;. To prove this, we proceed by induction. For 
A; = 1 we have (/o g)' = {f o g)g' , i.e., gi^i = g' and Lemma rA.2l implies gi^i G Sp-i{J). Assuming 
that the claim is true for k, we obtain 



X 



a—k 



k 



ifog)^'^ = Y,if^'^°9)9: 



k k 



if o g)^'+'^ = j;(/(^+^) o g)g^^,g' + ^^(/O) o g)g'^ 




= Y.(f^'^°9)hk+l 



where 



9i,k+i = g'l^k £ Sis^(^k+i)iJ) 

9j,k+i = 9j-i,k9' + g'j^k S Sji3^(^k+i)iJ), i = 2,3, . . . , A; 

gk+i,k+i = 9k,kg' s <S'(fe_|_i)^_(fc_,_i)(J) 



by Lemma lA.21 and the claim follows. Consequently, we have 




a(5-k 



□ 



We also state a simple corollary. 

Corollary A.4. Let f G So{I) and \f{x)\ ~ 1 for all xel. Then j G So{I). 

Finally, the inverse of a function that behaves like a symbol inherits this property. 

Lemma A. 5. Let a / 0, / G Sa{I) and suppose that \f{x)\ > |a;|", ^ |a;|"~^ for all x G I. 

Assume further that f'{x) ^ for all x e I. Then /"^ G Si (J) where J := f{I). 

a 

Proof. Since /' is nonzero on I, f~^ : J I with J := /(I) exists as a smooth function. By 
definition of the inverse and the assumptions on /, we have 

\y\ = \f{f-\y))\c^\f-\yr 

which implies < for all y & J. Furthermore, note that 

< ir'(2/)r-" < \y\^-' 

for all y G J by the assumption on /'. The claim now follows inductively based on this formula. □ 

Appendix B. Symbol behavior of solutions to Volterra equations 

In this section we discuss how the symbol behavior carries over to solutions of certain Volterra 
equations. We also allow the kernel to depend on an additional parameter A which is a situation we 
frequently encounter in this paper. It is then important to understand derivatives with respect to 
this parameter as well. We remark that the following result covers both, oscillatory and exponential 
behavior of the kernel. 

Proposition B.l. Let xq G M, Aq > and oj G C\{0}, Re(a;) > 0. Furthermore, assume that 
a(-,A) and b{-,X) are (possibly complex-valued) functions that satisfy 

\did'M^,x)\ < CkA^r'x-', \did',b{x,x)\ < Cfc,,(x)-"-'=A^-^ 

for all X > xq, a G (0, Aq) and G No where a > 1 and /3 > 0. Set 

K{x, y, A) := T e2'^"a(u, X)du e'^'^ybiy, A). 

Jx 

Then the equation 

POD 

<p{x,X)= K{x,y,X)[l + cp{y,X)]dy (B.l) 

J X 

has a unique solution <p{-,X) that satisfies 

|5ia,V(x,A)|<Cfe,,(x)-"+i-'=A^-^ 
for all X > Xq, X G (0, Aq) and k,i E Nq. 
Proof. We seek a solution A) of the Volterra equation 

^{x,X) = f{x,X)+ I K{x,y,X)(p{y,X)dy 

J X 

where 

POO 

f{x,X) := / K(x,y,X)dy. 

Jx 
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\{nny)\ = 



f'{r\y)) 



By the assumption on a and two integrations by parts, we obtain 



e^'^"a(u, X)du 



1 

2\uj\ 



e^^ya{y, A) - e^^^a(a;, A) - / e^^"9„a(n, \)du 



ry 

1+ / \dla{u, X)\du 

J X 



< p2Re(a;)j/ 



for all xq < X < y and A € (0, Aq) since Re(a;) > 0. This imphes 



\K{x,y,X)\ < 



-2Re{uj)y 



l%,A)|<|6(2/,A)| <(y)-"A^ 



for all xq < X < y and A G (0, Aq) by the assumption on b. Consequently, we obtain 



/ sup \K{x,y,X)\dy< / {y)-''dy<l 

Jxo {x:xo<x<y} J xq 



since q > 1, /3 > and a standard Volterra iteration yields the existence of a unique ip{-, A) satisfying 
Eq. dEH) and \ip{x, A)| < 1 for ah x > xq, X e (0, Aq). However, we have \ f{x, A)| < and 
thus, Eq. (|RT]1 imphes that in fact \^{x,X)\ < for ah X > Xq, X G (0, Aq). This settles 

the case k = £ = 0. 

In order to prove the derivative bounds note first that 



K{x, T] + x, X) 



Tj+X 



e2-"a(n, X)du e-2'^(^+^)6(r? + x, X) 
2^("+^)a(w + X, X)du e-2'^(''+^)6(?7 + x, A) 
e2'^"a(u + X, A)dn e'^'^nir] + x. A) 

JO 



and thus, 



£ k 



didlK{x,r^ + x,X) <Ck,eJ2Yl 



e'=ok'=o 

k 



a^'a(u + X, A)d'u e-2'^''a^-^'af-'='6(r/ + x, A) 



<Cfc,,A''-'5](:r)-'='(r? + x)-"- 



-fc+fc' 



fe'=0 



Now we proceed by induction in A: + ^. Let n £ N and assume that 



\did^^^{x,X)\ < Cfc,,,(x)-"+i-^''A'5-^' 
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for all k',n' e Nq with k' + 1' < n. We claim that this implies \d{dlip{x , \)\ < Cfc,^(j;)-"+i-'=A^-^ 
for all A;, f € No with k + 1 = n. Indeed, we have 

didl^{x,X) = J2Y. ii) ik) rdi-''dt''K{x,v + x,X)did'J[l + ^{v + x,\)]dv 
e=Qk'=o ^ ^ ^ ^ -^0 

/•oo 

/ di-''dt''K{x,v + x,x)d';d^J[i + ^{v + x,x)]dv 

, . «, Jo 



k'+e<n 

3^ ak 



+ / K{x,r] + x,X)d{d^(p{r] + x,X)dri 
Jo 

/•oo 

Oc((x)-"+^-'=A^-^) + / K{x, y, X)did'^^{y, X)dy 



and a Volterra iteration yields the claim. □ 

Furthermore, we discuss a variant which is useful when there is no decay in x at all. In this case, 
of course, oo cannot be chosen as a base point and one has to start from a finite xq. However, in 
certain situations it is possible to construct the solution on an interval which becomes infinite as 
an additional parameter tends to zero. 

Proposition B.2. Fix xq £ M, Xq > 0, a > —1, /3 > 7 > and assume that /3 — (a + 1)7 > 0. 

Let c be a real-valued function that satisfies c(A) > xq and |c(^)(A)| < CfX^'~^~^ for all X G (0, Aq), 
a S Nq. Furthermore, assume that the (possibly complex-valued) functions a{-,X), b{-,X) satisfy the 
bounds 

\did^a{x,x)\ < CkA^)-^^-\ \didlb{x,x)\ < CkA^r'"^"'' 

for all Xq < X < c(A), A G (0, Aq) and k,£ £ Nq. Set 

K{x,y,X):= j e-'^'^''a{u, X)du e'^'^Hiy, X) 
Jy 

for Xq < y < X < c(A) where u £ C\{0} and Re(a;) > 0. Then the equation 

/•X 

ip{x,X)= K{x,y,X)[l + ip{y,X)]dy 

JXQ 

has a unique solution f{-,X) that satisfies 

\did'M^,x)\ < CkAxr+'-'^x^-' 

for all Xq < X < c(A), A G (0, Aq) and k,£ £ Nq. 
Proof. As in the proof of Proposition IB. II we obtain the estimate 

\K{x,y,X)\<\b{y,X)\<{yrx'' 
for all Xq < y < X < c(A), A £ (0, Aq) and therefore, 

r-c(A) 

sup \K{x,y,X)\dy<X^-^^+'^^ <1. 



'xq {x:xo<y<x<c{X)} 

Consequently, a standard Volterra iteration yields the existence of (/?(•, A) with \if{x, A)| < 1 for all 
Xq < X < c(A), A £ (0, Aq) and we obtain 

Mx,X)\< r\K{x,y,X)\dy<{xr+'Xf' 

J xo 

in the relevant domain of x and A which settles the case k = i = 0. 



^The last condition follows from /3>7ifaG( — 1,0]. 
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Differentiating with respect to A we obtain 

d'Mx,X)=J2(l) r diK{x,y,\)di-''^{y,\)dy 

+ / K{x,y,X)diip{y,X)dy 

for any £ € N and the claim for A; = follows inductively by Volterra iteration. 

If fc > the situation is a bit more subtle since it is not completely obvious from the onset how to 
obtain decay in x for large x since the Volterra iteration starts at x = xq. The effect that provides 
decay can be illustrated by the following simple example: while Jq {y)~^ dy for N & N does not 
decay in x as a; — )• oo, we do have the estimate 



^^Hyr^'dy < Cn{x)- 



■N 



which follows immediately by means of one integration by parts and noting that y i— )• {y)~^ is 
monotonically increasing for y > N. In order to exploit this, one has to use some nice properties 
of the differentiated kernel 

Kix, y, A) := d,K{x, y, A) = e-^'^^^-y^a{x, A)6(y, A) 

which the original kernel K{x, y, A) does not enjoy. Consequently, the appropriate starting point is 
the equation 

dMx, A) = K{x,x, A)J 1 + ip{x, A)] + ^ K{x, y, A)[l + ip{y, X)]dy 



=0 



X — XQ 







k{x,x-ri,X)[l + Lp{x-r],X)]dr]. (B.2) 



Now we make the following two observations. 

(1) \d{dlk{x,XQ,X)\ < Ck/,N{x)-^ X^-^ for all xq < x < c(A), A € (0,Ao), k,i eNo and any 
N gN, 

(2) we have 



\did^,k{x,x- v,X)\ < Cfc,,A^-^e-2i^<^H^ - r/)"" 



-fc+fc' 

~ '// 

k'=0 

and thus, for any m G Nq, 



\did^,k{x, x-rj,X)\{x- r?)— dT? < C7fc,,,„(x)"-'=-'"A^-^ 

in the relevant domain of x, X and for all k,i £ Nq. 
Having this in mind we immediately infer from Eq. ()B.2p the estimate 

\dM^A)\ < {xrx^ + < (x)°A/3 

as claimed. Note that we have traded (x) for A^^ which is possible since /3 > 7 by assumption. It is 
now a straightforward matter to set up an induction based on Eq. ()B.2p that proves the remaining 
bounds for the higher derivatives. □ 
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Appendix C. Airy functions 

Apart from the results on Volterra equations with Airy kernels, the material in this appendix is 
mostly standard, see e.g [13], [TT] . 
Airy's equation 

y"{x) = xy{x) (C.l) 
is the simplest example of an equation with a turning point. It is well-known that the Airy function 

/>oo 

Ai(x) := M cos(\t^ + xt)dt, xGM (C.2) 
Jo 

is a solution of Eq. (jC.ip . The second standard solution of Airy's equation (jC.ip is traditionally 
denoted by Bi and defined by 

/•oo 

Bi{x) := 7 / [exp(-^t3 + xt) + sm{^t^ + xt)] dt 
^ Jo 

for X G M. We have the Wronskian relation l^(Ai, Bi) = -. By transforming the real integrals into 
contour integrals, one obtains the analytic continuations of Ai and Bi to the complex plane. The 
corresponding integral representations can then be used to study the asymptotics by the method 
of steepest descent. This is by now textbook material (see e.g., [TT] for a very nice account) and a 
classical result is 

Ai(x) = (47r)"5x~3e~§^^'^^[l + 0{x^^)] 

Bi(x) = 7r^5x^3ei^'^''^[l + 0{x^^)] 

as z — )• oo. Note, however, that this approach does not automatically yield the symbol behavior 
of the O-terms which is most important for our purposes. Consequently, we choose a different 
method based on the Liouville-Green transform. 

C.l. Asymptotics in the exponential regime. For x > 0, solutions to Eq. (jC.ip are expected 
to increase or decrease exponentially and thus, the idea is to transform Eq. (jC.ip into the form 

l/j" -1/; = Vip 

with a small right-hand side that can be treated as a perturbation. This is done by means of a 
Liouville-Green transform. Based on Section r2.lt we define the desired change of variables ^(x) by 

i^'(x)^ = X which immediately yields ^(x) = |x2. Clearly, : (0, oo) — )• (0,oo) is a diffeomorphism. 

Setting ^^{^{x)) = ^'(x)2y(x) = x^y{x), Eq. (IC.lj) transforms into 

<(e)-V'(0 = -irV(0 (C.3) 

and the right-hand side is indeed small for ^ large. Since we are interested in the asymptotics as 
X — >■ oo, we may restrict ourselves to ^ > 1. Note that the functions ^ i— t- e^^ are solutions to 
Eq. (jC.SP if the right-hand side is set to zero. Now we construct a fundamental system to Eq. (jC.SP 
by perturbing ^ i— )• e^^. 

Lemma C.l. Eq. (|C.3P has a fundamental system {'ip±} of the form 

V^±(0 = e±«[l + a±(0] 
where the functions a± are real-valued and satisfy the bounds 

for all and A; e Nq. 
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Proof. We start with ip-. According to Lemnia l4.1i it suffices to consider the Volterra equation 

poo pr) 

a_(0 = e'^'dv'e-'^{-l^)v-'[l + a_{r,)]drj. (C.4) 

Consequently, Proposition IB. II yields the existence of a_ satisfying the stated bounds. 

In order to construct the increasing solution -i/^+j we use the standard reduction ansatz, i.e., we 
set 



which is certainly well-defined for large enough c and > c and it is a solution to Eq. (lOSj) . If we 
define a+ by ip+iO = e^[l + CL+iO] it follows that 

a+(0 = 2e-2? [1 + a^iO] j^^ e^\\ + a.(r?)]d7? - 1 

where a_ satisfies |a^^(i^)| < Ck^,^^^^ for all ^ > 1 and G No (cf. Appendix [A]) . Consequently, 
we obtain 

o+ (0 = e-2« [1 + a„ (e)] (e2« - e^') + 2e-2« [1 + a„ (0] e'^'a- - 1 
= a_(e) - e2'=e-2«[l + a_(0] + 2[1 + a_(0]e-'« e2^a_(7?)dr?. 



Now we have 



e-2€ y e2''a_(7?)d7? = ia_(e)-ie2^e-2«a_(e)-ie-2« j e^'^d^{7])d7] = 0{C^) 

by noting that e'^'^d-{r]) is monotonically increasing for large tj and the 0-term behaves like a 
symbol by Appendix \^ This proves the claimed bounds for a+ and by solving an initial value 
problem, the solution can be smoothly extended to ^ > 1. □ 

Corollary C.2. For x > 0, the Airy functions Ai and Bi can be written as 

Ai(x) = (47r)"5(x)"3e~§'''^'"[l + a(x)] 

1 12 3/2 

Bi(x) = 7r"2 (x)~4e3^ [1 _!_ 5(2;)] 
where a and b are real-valued and satisfy 

for all X >0 and € Nq. Furthermore, Ai(x) and Bi(x) do not vanish for x > 0. 
Proof. Let ip± be as in Lemma I C. II and set 

y±{x) := x-3^±(^(x)) = x"3e^i'''^'[l + a±{^{x))]. 

By construction, {y±} is a fundamental system for Airy's equation (jC.ip and thus, there exist 
constants c± such that Ai(x) = c_y_(x) + c+y+(x). However, since Ai(x) approaches zero as 
rr — )• 00 whereas y+{x) grows unbounded, it is immediately clear that c+ = 0. Similarly, Bi must 
be of the form Bi(x) = d_7/_(x) + d+y+{x) for constants d±. This yields the stated representations 
for X > 1. However, the fundamental system {y±} can be smoothly extended to [0, 00) since Airy's 
equation is regular there. The fact that Ai(x) and Bi(x) do not vanish for x > is well-known and 
follows from the respective integral representations. □ 
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C.2. Asymptotics in the oscillatory regime. Due to the sign change of the right-hand side, 
solutions to Airy's equation (]C.ip oscillate on the negative real axis and therefore, their asymptotic 
behavior as x ^ — oo differs radically from Corollarv IC.2[ Indeed, the classical result [13j is 

Ai(-x) ± zBi(-x) ~ 7r"^e^*f a;~3e^*§^^''^ 

as X — )• oo. In order to estimate the error terms, it is convenient to switch from x to — x, i.e., we 
consider the Airy equation 

y"{x) = -xy(x) (C.5) 
on X > 0. Note that Ai(— x) and Bi(— x) are solutions to Eq. (jC.Sp . As before, we apply the 
Liouville-Green transform ^p(S,{x)) = xiy{(,{x)) with ^(x) = |x2 which yields 

i^"{o+m = -u-'m (C.6) 

and this time, the appropriate approximate solutions are ^ i— )• e^*^. 
Lemma C.3. There exists a fundamental system for Eq. ()C.6p of the form 

V'±(0 = e±*«[l + a±(0] 

where the functions a± satisfy 
for all and A; e Nq. 

Proof. The existence of ip- with the stated properties follows from Proposition IB. II and "0+ is just 
the complex conjugate of tp^. □ 

As a consequence, we immediately obtain the following representation of AiitiBi in the oscillatory 
regime. 

Corollary C.4. For x > 0, the Airy functions Ai it iBi have the representation 
Ai(-x) ±iBi(-x) = 7r-^e±^f (x)-ieT*§^'^'[l + a±(x)] 

where the a± satisfy 

|af^(x)| <Cfc(x)-i-'^ 

for all x >0 and k gNq. 

C.3. Volterra equations with kernels composed of Airy functions. For equations with 
turning points, the Airy functions are of fundamental importance. Consequently, they frequently 
appear in the kernels of Volterra equations in perturbation theory. The following result is the 
analog to Proposition IB. II in the exponential case uj = 1. 

Proposition C.5. Fix Aq > and suppose a{-,X), b{-,X) are (possibly complex-valued) functions 
that satisfy 

\did'MxA)\ < Cfc,,(x)-'^A-^ \did',b{x,X)\ < Cfc,,(x)-"-'=A'^-^ 
for all X >0, Xe (0, Aq) and k,ie'No where a > |, /3 > 0. Set 

rv 

K{x,y,X):= / Ai{uy'^a{u, X)du Ai{yfb{y, X). 



Then the equation 

poo 

ip{x,X)= K{x,y,X)[l + ip{y,X)]dy 

J X 

has a unique solution (p{-,X) that satisfies 

|ai5^<^(x,A)| <C,,,(x)-"+5-^A'^-^ 
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for all X >0, X e (0, Aq) and k,i e Nq. 

Remark C.6. We emphasize two important differences between tfie exponential and tlie Airy case. 
First, the required decay in x of the function b{x, A) is weaker. In the Airy case one only needs 
a > ^ whereas in the exponential case one has to require a > 1 (see Proposition IB. ip . Second, the 
loss of decay in x of the solution if{x, A) compared to b{x, A) is weaker. In the Airy case one only 
loses (x) 2 as opposed to (x) in the exponential case. 

Proof of Proposition IC.51 The strategy is to reduce the problem to Proposition IB. II Define rj : 
M — 7- R by 7]{^) := Then r/ is a diffeomorphism and we have r]{^) = (|)3^3 for ^ > 1. 

Furthermore, we note that t/'(^) ~ (0^^ and \r]^^'^{C)\ < (^^=(0^"^ foi' all ^ G R and A; G No- We 
intend to study the equation 

ip{x,X) = / K{x,y,X)[l + ^iy,X)]dy= / ry(0, A)[l + <^(r/(e'), A)]r?'(0'ie'. 

Hence, writing (^(^, A) = f{r]{^), A), we obtain 

/•oo 

^ic,x) = j^ K{c,e,m+^ic',mc' (C.7) 

where K{^,(^',X) := K (rj^^) , r](^') , X)rj' (^') . By assumption, the kernel K is of the form 

pin 

Kit e', A) = / Aiiuy^aiu, X)du Aiir^iOfbir^iO, >^Wi^') 

= Mivinr'airjiaAWiadC" Ai(r?(e'))'K^(^'),A)7?'(e')- 

According to Corollary IC. 21 we have the representation 

Ai(r/(e))±2 ^ (4nr\r^(^))TU'f'^[l + a±itX)] 
where \d{d^a±{t X)\ < Cfc.XO^^^^'A-^ for alU > 0, A G (0,Ao) and k,£ G Nq. Writing 

a(e, A) = a(77(0, ^WiOiviO) ^ [1 + 

A) = biviO, >^WiO{m)'Hl + a+iO], 

we obtain the bounds 

\didld{tX)\ < CkAO'"^-', \didfb{C,X)\ < Cfc,,(0-^°-^'A/^-^ 
for all ^ > 0, A G (0, Aq) and k,£ £ Nq, see Lemma |A.3| and K can be written as 

me,X) = l^^ e'^"d{e',X)dC" e-^«'6(e',A). 

Consequently, Proposition IB. II applies to Eq. ()C.7p and we obtain the bounds 

|5l9f(^(e,A)| <Cfc,,(0-i"+^-'A^-^ 
which, via Lemma lA. 31 translate into the stated bounds on ip. □ 
Next, we consider the oscillatory case. 
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Proposition C.7. Let Aq > and suppose a(-,X), b{-,X) are (possibly complex-valued) functions 
that satisfy 

Idid'Mx, A)| < Cfc,,(x)-'=A-^ \did',b{x, A)| < Ck^eixr^'-'X"-' 
for all X >0, X £ (0, Aq) and k,i eNo where a > ^, /3 >0. Set 



K{x, y, A) := / [Ai(-u) ± iBi{-u)]-^a{u, X)du [Ai(-y) ± iBi{-y)fb{y, A). 

Jx 

Then the equation 

POO 

ip{x,X)= K{x,y,X)[l + ip{y,X)]dy 

J X 

has a unique solution (p{-,X) that satisfies 

|aia^<^(x,A)| <C7fc,,(x)-"+5-'=A'^-^ 
for all X >0, X £ (0, Aq) and k,i £ Nq. 

Proof. Based on the representation of Ai it iBi from Corollary IC.41 the problem can be reduced to 

POO 

(^(e,A) = y^ i^(e,e',A)[l + (^(e',A)], ^>0 

where 



k{^, e', A) = 1 e^'^^"a{e, X)de e^^ni^, A) 



and d, b are like in the proof of Proposition IC.5I Thus, as before, the claim follows from 
Proposition IB.li □ 

Finally, we state the "Airy version" of Proposition IB.21 

Proposition C.8. Fix xq G M, Aq > 0, a > /3 > §7 > and assume that (3 — {a + > 0. 
Let c be a real-valued function that satisfies c(A) > xq and |c*^^)(A)| < CfX^'~^^^ for all X G (0, Aq), 
H. € No- Furthermore, assume that the (possibly complex-valued) functions a{-,X), b(-,X) satisfy the 
bounds 

idid'Mx, A)i < Cfc,,(x)-'=A-^ \did',b{x, A)i < CkA^r-'x^-' 

for all xq < x < c(A), A G (0, Aq) and k,£ £ Nq. Set 

ry 

K{x,y,X):= / Biiuy^aiu, X)du Bi{yfb{y, X) 



for xq < y < X < c(A). Then the equation 

(p{x,X)= K{x,y,X)[l + <p{y,X)]dy 



XQ 



has a unique solution (p{-,X) that satisfies 

|aia,V(3^,A)| <Cfc,,(x)"+^'=A^-^ 
for all xq < X < c(A), A G (0, Aq) and k,£ £ Nq. 

Proof. Analogous to the proof of Proposition [C31 the Volterra equation for if we want to solve can 
be rewritten as 
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where t]{^) := {^)h{0~K ^{t>^) ■= Co := ry-^^^o) and 



e, A) = f e-2«"a(e", A)^^" e2«'6(e', A). 



The functions a and b are given exphcitly in the proof of Proposition IC.5I in terms of a and b and 
they satisfy the bounds 

\didl~ai^,X)\ < Cfe,,(0-'=A-^ \didfbitX)\ < Cfc,,(0^"-^"'A^-^ 

in the relevant domain. Furthermore, we have the bounds \dir]-\c{X))\ < CeX'^^'^ for all A G 
(0,Ao). Consequently, Proposition IB.2I vields the existence of satisfying the bounds 

|9i9f(^(^,A)| <Cfe,,(Oi"+^-'A'^-^ 
and the claim follows. □ 



Appendix D. Modified Bessel functions 

For the convenience of the reader we provide the essential details of the asymptotic theory for the 
modified Bessel functions. The following results are mostly standard, see Chapter 7, §8, Chapter 
10, § 7, and Chapter 11, § 10 of [13], as well as [1], [9], [19]. However, as before in Appendix [Cl 
our focus is a bit different and we derive the necessary results based on the theory developed in 
Appendices \M M and O 

Consider the modified Bessel equation 

w^y"{w) + wy'{w) + (z/^ - w'^)y{w) = (D.l) 

where in general w,!^ G C, however, for our purposes it suffices to consider w > 0, G M. In 
particular, we are interested in the behavior for large. Solutions of Eq. (jP.ip can be obtained 
by a standard Frobenius ansatz around the regular singular point w = which yields the modified 
Bessel function li^, given by 

where the function b satisfies the bounds \dld^b{w , 1^)1 < Ck/w^^^^'^^'^'^^ for, say, all < to < 2, 
< 1 and k,£ £ Nq. In order to find the asymptotics of Iiy{w) for w — oo, one needs some kind 
of global representation of Iiy which is provided by the well-known integral formula 

h.{w) = (1 - tT-'^e^'dt. (D.3) 

■K2T{\ + iu) J -I 

The classical result for the asymptotics of li^ is 

Ii^{w) = (2™)-^ [e"'(l + 0{w-^)) - ie^"e-"'(l + 0{w-^))\ (D.4) 

as — )■ oo. 



D.l. Asymptotics for small v. For our purposes it is important to be able to construct solutions 
to the modified Bessel equation via Volterra iterations. This procedure has the advantage that it 
automatically yields derivative bounds for the error terms. Furthermore, it may be considered as 
a warm-up exercise for the large v asymptotics studied below and also, for the more complicated 
problems we have to deal with in this paper. 
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As always, it is convenient to remove the first derivative in Eq. (jP.ip by defining a new dependent 
variable B{w) := w^y^w) which consequently satisfies 



At this point it should be noted that are solutions to Eq. (jP.SP and from Eq. ()D.2p we 

immediately obtain 

W{V-I-iu, V-Ii.) = f sinh(7ri^) 

which implies that \f-Iiu} is a fundamental system for Eq. ()D.5P provided that G M\{0}. 

In particular, this shows that there is no > such that I-i^iwo) = liuiwo) = and, since 
lii, = I-iu, we conclude that Iip{w) does not vanish for any w > and i/ G R (the case = 
follows directly from Eq. (|D.3p ). 



D.1.1. The exponential regime. We proceed by constructing a fundamental system for Eq. (jP.SP 
with known behavior as — )• oo. This is the exponential region of the equation and we expect to 
obtain an exponentially decaying and an exponentially increasing solution. Matching this system 
to ^/'■Iiu yields bounds for the derivatives of the O-terms in Eq. ()D.4p . In order to construct the 
desired fundamental system, we rewrite Eq. (ID.5j) as 

B"{w)-B{w) = ^B{w) 



and treat the right-hand side perturbatively. 

Lemma D.l. Fix t^o > sufficiently large and z^o > arbitrary. Then there exists a fundamental 
system {Bi{- v) , B2{- ., v)} for Eq. ()D.5P of the form 

Bi{w,v) = e-'"[l + hi{w,v)] 
B2{w,v) = e'"[l + b2{w,v)] 

where the functions bj{-,i'), j = 1,2, are real-valued and satisfy 

\didtbj{w,iy)\<Ck,iW''-' 

for all w > wq, V G [— I'Oi '^o] ^f^^^ k,i €z Nq. 

Proof. We start with the decaying solution Bi. According to Lemma |4.H we have to construct a 
solution bi{-,i') to the equation 

foo fv / w2 _|_ 1 \ 

b^(w, v)= / e^'^du ^ [1 + u)]dv 



w J w 



for w > wq and thus, Proposition IB. II yields the desired result. If wq is sufficiently large, we 
obviously have Bi{w,v) > for all w > wq and we obtain B2{-,u) by the standard reduction 
ansatz, i.e.. 



B2{w,u) := 2Bi{w,u) 



/ Bi{v,iy)-^dv + 1 

J wo 



cf. the proof of Lemma IC.ll □ 
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D.1.2. The oscillatory regime. Next, we deal with the oscillatory regime. The behavior of solutions 
in this regime is completely described by Eq. ()D.2p and there is no need to employ Volterra iterations 
since everything follows by Frobenius' method. As a consequence, we obtain the following. 

Corollary D.2. Fix Wo,i'q > 0. Then there exists a fundamental system {B±{-,i')} for Eq. ()D.5P 
of the form 

B±{w, v) = w^^'''[l + b±{w, v)] 
where the functions h±{-,v) satisfy the hounds 

for all < w < wq, < < z^q cind k,£ £ Nq- 

D.1.3. Matching of the fundamental systems and the classical Bessel asymptotics. In order to be 
able to glue together the fundamental systems {B±{-,i>)} and {Bj{-,iy) : j = 1,2} obtained in 
Lemma ID. II and Corollarv ID. 21 we have to calculate the Wronskians W{B±{-,i'),Bj{-,i')), or, 
equivalently, the precise asymptotics of B±{w,i') as w — )• oo. However, without any global infor- 
mation, one can only obtain a trivial result here, namely that the Wronskian is a smooth function 
of u. Thus, at this point we have to make use of the global information provided by the classical 
asymptotics Eq. ()D.4p that follow from the integral representation ()D.3p . Note that the following 
result yields the desired symbol character of the 0-terms in Eq. ()D.4p . 

Lemma D.3. Let wq > ^ he sufficiently large and uq > he arhitrary. The Bessel functions 
B-i-^-ju) from Corollary \D. 2 have the representation 

B±iw, u) = 2±^'^r(l ± iu)^I±i^{w) 

= ^-52-5±*'^r(l ± iv) [e^'il + h2{w, ±y)] - ie^'^^e-'^il + bi{w, ±iy)]] 

where the error terms 6j(-,ibz^), j = 1,2, are real-valued and satisfy the bounds 

\didtbj{w,±i^)\<C,,,w-^-'' 

for all w > Wq, < < z^o ^.i^d k,i £ Nq. In particular, we have 

W{B±{-, u),Bi{-, ly)) = -7r"^25±»'^r(l ± iu) 

and \B-i-{w, z^)| > for all w > and < < vq. 

Proof. According to Eq. ()D.2p . we have the asymptotics 

2=F*^ 



wI±Uw) ~ " . W2^''' {w ^ 0+) 

1 (1 ± lU) 

and thus, by comparison with the asymptotics of i^), we obtain 

B±{w, v) = 2^'''T{1 ± iv)^I±,^{w). 

This implies \B±[w,v)\ > for all u; > 0. For notational convenience we set B±{w,v) : = 
(27r) 2y^/_i_j^(ti;). From Eq. ()D.4p we have the asymptotics 

B±{w, v) = e'"\\\ 0{w-^)\ - ie±^"e-"'[l + 0{w~^)\ 

as u; — )• (X). Since i?±(-, is a solution to Eq. (ID.5P and i^) ■ j = 1, 2} from Lemma fD. II is a 

fundamental system for that equation, there must exist connection coefficients such that 

B±{w, u) = ^±^i{u)Bi{w, u) + i3±^2{i^)B2{w, v). 

Dividing this equation by and letting w — )• oo, we obtain (3±^2i'^) = 1- This implies 

lmB±{w,i') = Im/3±^i{i')Bi{w, v) 
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and multiplication by e"' yields lm.f]±^i{iy) = —e^^'" in the limit w — )• oo. Since B±{-,v) = 
{2'k)-^22^^''T{1 ±iv)B±{-,u), we infer 

1 W{B±{-,u),B,{;v)) _ ~ . ^ _ 

(2^)-|2±-r(l±ii.) W{B^{.,u),B2{;u)) ^^''^ ' 

and the claim follows from W{Bi{-^ v),B2{-, v)) =2. □ 

D.2. Asymptotics for large v. For the construction of fundamental systems of Eq. (jP.Sp in the 

case of large a more careful treatment of the turning point is necessary. The character of the 
solutions is dictated by the sign of (i^^ + | — i.e., if i^^ + | — it;^ > solutions tend to oscillate 
whereas for i^^ + 4 — ij;^ < we have exponentially increasing or decreasing behavior. The turning 



point wt which separates these two regimes is wt{y) = \j + \ and it moves with u and behaves 

like wt{v) ~ V for — )• oo. Consequently, it is reasonable to base the analysis of Bessel's equation 
on Airy functions. However, it turns out to be crucial for the asymptotic analysis to fix the turning 
point first by introducing the rescaled variable x := ^, i.e., we set ip{x) := B{vx) and Eq. (1D.5P 
transforms into 

^"(^) _ ^2^1 _ J^)^(^) = -^^(x). (D.6) 

on X > 0. Heuristically speaking, we expect the right-hand side to be negligible provided that v 
is sufficiently large. Consequently, we ignore the right-hand side altogether for the moment and 
observe that the turning point for the "remaining equation" is now fixed at x = 1. 

D.2.1. Reduction to a perturbed Airy equation. The idea is now to define new variables such that 
Eq. (|D.6P turns into a perturbed Airy equation. This is accomplished by means of a Liouville-Green 
transform. According to Eq. (j2.2p we require 

C'ixfa^) = 1-^ (D.7) 

which yields 



= sign(x — 1) 



1 



dy 



fD. 



by separation of variables, cf. Lemma [3.11 Note that C is chosen in such a way that the turning point 
X = 1 is mapped to 0. By a Taylor expansion around x = 1 it is readily seen that C, : (0, oo) — )• M 
is smooth (in fact, real-analytic) and obviously surjective. Furthermore, we have C' > and 
this shows that C maps (0, oo) to M diffeomorphically and thus, C defines an admissible change of 
variables as desired. Consequently, Eq. ()D.6P transforms into 

f{c) = iy\m +y2{c)m^ (d.9) 

Airy pert, 
for C e M with ^(((x)) := C {x)^2^l,{x) and, cf. Eq. ^1?^ . 



3 C"{xf 1 C'"(x) 1 5 C(^)x2(4 + x2) 

V2{ax)) = -7^773T + o777-T?^ " ^^'?.^,f^^'>. = 71?773^ + 



4C'(x)4 2C'{xY 4x2C'(x)2 16C(x)2 4(l-x2)3 

where the last expression follows straightforwardly from Eq. ()D.7p . Note that the apparent singu- 
larity at X = 1 cancels by a Taylor expansion and the fact that C(l) = 0. Consequently, we have 
V2 G C°°(M) and it is evident from the above expression and Appendix IA1 that V2 satisfies the 

bounds \ V^^\0\ < Ck{C)~'^~'' for all /c G No and aU (gR. 
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D.2.2. The exponential regime. The Airy equation 

<P"iO = y\m 

2 2 

has the fundamental system {Ai(i/3-),Bi(i/3-)} where Ai, Bi are the standard Airy functions, see 
Appendix [O This immediately allows us to construct a fundamental system for Eq. (ID. Oh in the 
exponential regime. 

Lemma D.4. For C > 0, 1 there exists a fundamental system {</>i(-, i/), (/)2(-, i^)} of Eq. ()D.9P 

of the form 

<Ai(C,J^)=Ai(i/§C)[l + ^-'ai(C,^)] 
</.2(C,z^) = Bi(z.§C)[l + i^"'«2(C,'^)] 

where the functions aj, j = 1,2, are smooth, real-valued and \aj{C,i')\ < 1 m the above range of Q 
and V . Furthermore, aj satisfy the bounds 

|9^9fa,(C,^^)|<Cfc,,(C)-i-V-^ C>1 

as well as 

\diajiO,i^)\ < C,v-\ \didi^aj{Q,v)\ < C^J-' 
for all G No and > 1. 

Proof. We start with cpi and according to Lemma l4.ll we have to construct a solution oi of 

J Ai{iy3r]'y^dr]'Ai{u3i]fV2{r])[l + v'^ai{ri, u)]dr] 



I I L Ai{u)-'^duAi{yfV2{u--^y)[l + u~^ai{u--iy,u)]dy 

Ji/'SC J u^C 



for C ^ 0. Thus, the function ai{x,v) := u~^ai{i'~3 x) satisfies the Volterra equation 

ai{x,u)=j K{x,y,u)[l + ai{y,u)]dy 

J X 

where 

ry 4 2 

K{x,y,v) := I A:\{uy'^du A\{y)'^v~^V2{i>~^y). 



In view of Proposition IC.5I we set X := u ^ and, since we are interested in large v, we assume 
A G (0, 1). If X > A^3 > 1, we have x = (x) and hence, (Aax)"^ = A^sx^^ which yields the bounds 

\did^,\IV2ixlx)\ < Cfe,,(x)-2-^•A-^ fc,^ E No 

2 

for all x > A~3 and A G (0, 1). Consequently, Proposition IC.5I yields the existence of ai with the 
bounds 

\did%ix,\-')\<C\A^)-'l-''X-' 
and, via Lemma \A.3\ these translate into 

\did%{x,u)\<CkA^r'^'''^'' 

2 2 2 3 

for all X > z^3 , zy > 1 and k,£ G Nq. We obtain |ai(C,z^)| = |i^ai(z^3 z^)] < v{v^C)~2 and 

\didlar{Cv)\<CkAC)-'^-''y-' 

41 



2 2 

for all > 1, z/ > 1 and k,i £ No as claimed. In the case < a; < A~3 (which implies (x) < A"^) 
we have 

2 

2 r^~^ 

ai(x,A-i) = ai(A-3,A-i)+ / K{x,y, X-^)[l + di(y, X-^)]dy 

J X 

and from the above we obtain the bounds |(9^ai(A^3 , A^-*^)! < C^X^^^ for i G Nq. Furthermore, we 
have 

|aia^ty2(Aix)| < c,,,(x)-^'At-^ 

2 

for all < X < A^3, A G (0,1) and fc,^ G No and, similar to Proposition IC.Sl (cf. also Proposi- 
tion IC.Sp , we infer inductively by Volterra iterations that 

|5ia^ai(x,A-i)<Cfc,,(x)-'=Ai-^ 

which implies the claimed bounds for ai(0, v) and 5^ai(0, v). 

For the function (/)2 we apply the standard reduction ansatz, i.e., we use the fact that (/>i(Ci v) > ^ 
for all > provided that v is large enough (hence the condition ^ 1) and set 



By recalling that VF(Ai, Bi) = ^ it is readily verified that is indeed of the stated form. □ 

D.2.3. The oscillatory regime. For < we cannot apply the same construction as in Lemma lD.41 
due to the zeros of Ai and Bi on the negative real axis. Consequently, we construct a complex 
fundamental system based on Ai it iBi. 

Lemma D.5. For C ^ 0, z/ > 1 there exists a fundamental system {(/)±(-,i^)} of Eq. (1D.9|) of the 
form 

MC, ^) = [Ai(i^K) ± iBi(z.§C)][l + ^-^a±{C, ^)] 
where the functions a±{-,v) are smooth and \a±{C,,v)\ < 1 in the above range of Q and u. Further- 
more, a± satisfy the hounds 

\didla^{C,u)\ < Cfc,,(C)-^"'^"', C < -1 

as well as 

|5^a±(0, 1^)1 < C,u-', |a^aca±(0, v)\ < C^^.v^-' 
for all I' > 1 and k,i G No- 

Proof. After switching from to — C, the proof becomes completely analogous to the proof of 
Lemma ID.4I The only difference being that one has to apply Proposition IC.7I instead of Propo- 
sition IC.Sl Furthermore, the reduction ansatz is not needed: once (p- is constructed, one sets 
0+ := JZ. □ 

D.2.4. Matching of the fundamental systems. Finally, we match the fundamental systems {4>j{-, v) ■ 
j = 1,2} and {(p±{-,i')} from Lemmas ID. 41 and ID. 51 This is done by evaluating their Wronskians 
at C = 0. Unlike in the small z^ case, we do not use any global information on the Bessel functions. 
This is not necessary here since we have a small parameter z^~^ and the Airy functions are defined 
globally. 

Lemma D.6. For C ^ 0, the functions (j)±{-, v) from Lemma \D.5\ have the representation 

</'±(C, i^) = [1 + ^~'^a±,i{'')]MC, ^) ± i[l + i^~'a±,2(z^)]<^2(C, ^) 
where (pj, j = 1, 2, are from Lemma D.4\ and the coefficients a±j satisfy the bounds 

\dia±j{u)\ < C,u-' 
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for a// > 1 and £ G Nq, j = 1, 2. 



Proof. The ansatz 
implies 

"^•'^''^ VF(c/>i(-,i.),</.2(-,i^))' "^•'^''^ Ty(02(-,i^),./'i(-,z^))' 
Note that the Wronskians are functions of v only since Eq. ()D.6P does not have a first derivative. 
Setting ^1 := Ai, A2 := Bi and A± := Ai it zBi, we obtain, for j = 1, 2, 

(^,(0, z.) = ^,(0)[1 + 0{v-^)] </.;.(0, u) = i.i4.(0)[l + 0(1.-1)] 

<A±(0, z.) = ^±(0)[1 + Oc{v-')] </.'±(0, 1.) = z.i^^(0)[l + Oc{v-^)] 

where all 0-terms behave like symbols. This yields 

a+ \(v] = 5 = 1 + Ov[v ) 

vzW{A^,A2)\\ + 0(v-^)\ 

provided that v is sufficiently large, and, analogously, 

a±,2(z^) = 2 = ±z + C'c(i. ) 

v-m{A2,A{)\\^0(v-^)\ 

where the Oc~terms behave like symbols. □ 
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